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Abstract—The selective properties of the ion-permeable membranes “ Permaplex A-10” 
and “ Permaplex C-10 ” in concentrated sodium-chloride solutions have been studied. Transport 
numbers have been obtained as a function of the mean concentration and concentration ratio 
existing outside the membranes. For a given mean concentration, both membranes show higher 
selectivity for high concentration ratios. Membrane transport numbers approach 0-6 at saturation. 


It has been shown that a three compartment cell, operating under special conditions, may 
demonstrate the true selectivity of the membranes. Back diffusion of electrode products was 
limited by high flow rates and short running times. Current efficiencies fall below those predicted 
from transport numbers, largely due to electro-osmosis, and are less than 20% at saturation. 


Résumé— L‘auteur étudie les proriétés sélectives des membranes semi-perméables “* Permaplex 
A-10 " et * Permaplex C-10 " dans des solutions concentrées de chlorure de sodium. Il a obtenu 
les nombres de transfert en fonction de la concentration moyenne et du rapport des concentra- 
tions a l'extérieur des membranes. Pour une concentration moyenne donnée, les deux membranes 
montrent une sélectivité plus grande pour des rapports de concentration élevés. A la saturation, les 
nombres de transfert s'approchent de 0-6. 


L’auteur montre qu’aveec une cellule 4 trois compartiments travaillant dans des conditions 
déterminées, on peut mettre en évidence la séléctivité réelle des membranes. Grace a des vitesses 
d’écoulement rapides et des durées suffisamment courtes, on peut alors limiter la diffusion en 
retard des produits aux électrodes. Les rendements en courant descendent au-dessous des 
valeurs prévues & partir des nombres de transport, en grande partie & cause de I'éctro-osmose. A 
la saturation, ils sont réduits de 20%. 


had to reproduce conditions analogous to those 
within such a cell. 

The general properties of ion-permeable 
membranes and the theory of their operation 
have been widely discussed and will not be 
repeated here. So_uner [15] has given many 
references to earlier work on membrane prepara- 


1. INTRODUCTION 


Tue present research was undertaken to 
determine the feasibility of applying ion-exchange 
membranes as separating media in concentrated 
brine solutions. The “ Permutit ” membranes 
used are typical of the few types now produced. 
They combine a high degree of selectivity with 


minimum electrical resistance and will thus give 
maximum separation with minimum power 
usage. The experiments were therefore under- 
taken for the final purpose of making design 
estimates for a large-scale, multicompartment 
separating unit, and the apparatus necessarily 
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tion and theory. 

Commercial membranes of the type used in 
this research have been described by Jupa 
and McRae [9], Kressman [11], WINGER et al. 
[20] and “Anon.” [2]. Technical data on the 
the membranes “ Permaplex C-10” and 


“ Permaplex A-10” have been given by the 
Permutit Co. Ltd. [13]. These commercial 
membranes have been used both for analytical 
and laboratory purposes (WILLARD and FINLEY 
[19], Friscu [3, 4, 5], ANpERsoN and Wy.am [1], 
Suarr et al. [14], and Kunin [12]) and for desalting 
in multicompartment cells (WincER et al. [21] 
Watrers, Weiser and Marek [17], Wrecners 
and Van Hoek [18], and Ionics Incorporated 
[8}). 


2. EXPERIMENTAL 
(a) Membrane transport numbers 


The transport numbers of the selected ion in 
Permaplex C-10"" and “ Permaplex <A-10” 
were computed from concentration potentials 
existing across the membranes when separating 
solutions of different activity. 

Consideration of molar free-energy changes 
involved in such a cell leads to the expression 


E = RT /F (2t, — 1) Ina,*/a,* (1) 
where 


E = membrane potential 
F = Faraday constant 
t, = transport number of selected ion 
4, = mean ionic activities in solutions 
1 and 2 


The theoretical maximum value of the concen- 
tration potential is given by 


E, = \na,*/a,*. (2) 


Thus 


t, =(E, + E)/2E,. (3) 


Equation (3) gives an average transport 
number in the range a, + ag. If the activities 
of the solutions surrounding the membrane are 
almost equal, ¢, may be assumed linear in In a, ay 
and will then be a function of the geometric 
mean of a, and a,. Other investigators have 
used various parameters to express concentration. 
SoLLNER [15] used the upper concentrations with 
2/1 ratio. Wincrr, Bopamer and Kunin 
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[20] used the arithmetic mean concentration with 
the same ratio. Jupa et al. [10] have used a 
ratio of 2/1 and plotted transport numbers 
against geometric mean concentration. For 
this work an arithmetic mean was arbitrarily 
chosen, since such a plot will reveal whether 
the concentration gradient through the membrane 
is linear for various external-concentration 
ratios. 


(i) Procedure. Membrane potentials were 
measured in a flow cell of the type given below : 


The cell was made in two halves from welded polyethylene, 
and assembled with the membrane between, acting as 
its own gasket. The solutions in use were fed at a rate 
of 30 ml/min to give a complete volume change in under 
one minute. Potential measurements were made using a 
Leeds and Northrup type 7552 potentiometer (K-2) with a 
galvanometer type 2430 of the same make. When two 
consecutive readings taken at intervals of one minute 
agreed within 0-05 mV, the potentiometer reading was 
taken as the concentration potential under the conditions 
applied. The symmetry of the electrodes was checked 
before and after every potential measurement. Readings 
were made at a temperature of 25-0C, with arithmetic 
mean concentration values ranging from 0-06 molar to 


2.88 molar. The molarity ratios studied were 2, 4, 8, 
and 20 to 1. 
(ii) Results. The transport numbers for the 


sodium ion in “ Permaplex C-10” and chloride 
ion in “* Permaplex A-10” are plotted in Fig. 1 
and 2 as a function of the arithmetic mean con- 
centration. Figs. 3 and 4 show transport numbers 
as a function of concentration ratio and indicate 
how the curves for a 1/1 ratio were obtained 
(see Figs. 1 and 2). 

Activity data were taken from Harnep and 
Owen [7] and the data for the saturation line 
in Figs. 1 and 2 from Handbook of Chemistry 
and Physics (6). 


(b) Three-compartment membrane cell 


The three-compartment cell is the simplest unit 
demonstrating the properties of both anion and 
cation selective membranes. Friscu [8] used 
the cell for electrolytic precipitation of uranium. 
ANDERSON and Wy tam [1] demineralized sugar 


Solution Sat. KCl Hy/Hg, Cl, 
M, Agar 
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solutions in a_ three-compartment cell with 


Permaplex membranes, while WILLARD and  petween the membranes of a three compartment 


FinLey [19] used a similar cell for separation of ¢e]] to obtain decontamination of radioactive 
uranium tetrafluoride. Watrers, Weiser and wastes. 


Marek [17] proposed placing ion-exchange resin 


1-0 


° 
© 


4 8 20 
Molarity ratio M, /™, 


Chloride -ion transport number 
@ 


° 


° 1 2 3 4 , 5 Fic. 4. The effect of concentration ratio on anion trans- 
Arithmetic mean molarity NoC! “ne past theough Formagion A-10. 
; (1) 0-5 Molar (4) 2-0 Molar 
Fic. 2. Transport number of chloride ion in Permaplex (2) 1-0 Molar (5) 2-5 Molar 
A-10. (3) 1-5 Molar (6) 3-0 Molar 


_ | 
99 


G. G. Mapewick 


Mass transfer occurs through ion-permeable 
membranes by four mechanisms : 
(i) normal electrodialysis, 
(ii) ionic diffusion, 
(iii) osmosis, 
(iv) electro-osmosis. 


Transfer by mechanism (i) may be predicted 
from the membrane selectivity : 

Consider a three-compartment cell containing 
in all parts a solution of the monovalent salt MX. 
(The various compartments need not contain 
solutions of equal strength). The centre compart- 
ment is bounded by a cation membrane on the 
cathode side and an anion-permeable membrane 
on the other. A potential is applied to the system 
through suitable electrodes in the outer compart- 
ments. 

In a given time interval let M* cations pass 
through the cation membrane towards the 
cathode. At the same time 2~ ions will pass in 
the opposite direction due to the imperfect 
nature of the membrane. Also in this time 
interval, let X~ anions pass towards the anode 
and m* ions in the opposite direction. 

The efficiency of the cation membrane 


net ions passed to catholyte 

total ions transferred 
M* 

which by inspection 

— ty 

+ 

= 


Ne 


where 
t.” = transport number of cation in cation 
selective membrane 


{> = transport number of anion in cation 
selective membrane. 


Similarly, the efficiency of the anion membrane 
= — 1. 


The cell efficiency is the arithmetic mean of 
Ne and 


This efficiency, which ignores mechanisms (ii), 
(iii), and (iv), will in general be different from that 
of a real membrane under operating conditions. 
Diffusion and osmosis have an additive effect 
equivalent to a net transfer of solute from con- 
centrated to dilute solution adjacent to the 
membrane. In a three-compartment cell, electro- 
endosmosis and electro-osmosis cause solvent 
transfer towards the electrodes at a rate dependent 
on ambient concentration, applied potential, 
and the distribution of charged groups in the 
membrane. 


(I) Apparatus. The three-compartment cell used is 
shown in Fig. 5. The electrode chambers were machined 
from methacrylate plastic, with cemented joints. The 
centre compartment spacer was cut from polythene 
3mm thick and was baffled to give six passes across 
the membrane face, and excellent mixing throughout the 
compartment. (Flow patterns in the compartment were 
examined by dye injection before a final baffle layout 
was made.) The cell had an overall diameter of 20 cm 
and a free membrance area of 100 sq. em. Liquid holdup 
was kept low and was about 25 ml for the centre and 
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100 ml for the electrode compartments. By placing a 
perforated support plate outside the membranes and a 
positive pressure in the centre compartment, disturbance 
of the membranes during operation was avoided. Power 
was supplied from a rotary d.c. generator to graphite 
electrodes of face area 80 sq. cm. Conventional rheostat 
control with voltmeter, ammeter, and a copper coulometer 
was provided. 

A flow diagram of the arrangement is shown in Fig. 6. 
Originally it was intended to run the electrode feeds at 
high rates and to follow the change in concentration 
of the recycled centre stream with a flow-type conductivity 
cell (Philips-type G.M. 4247, bridge-type G.M. 4249). 
This system was unsuitable since during a running time 
of 30-60 min the back diffusion of electrode products 
(ie. HgO* and OH™~ ions) caused pH changes which, 
although not significant in terms of ionic transfer, had a 


6 6 


Fic. 6. Flow diagram of three-compartment cell. 
(1) Cell (5) Flow meter 
(2) Electrode feed (6) Solution storage 
(3) Head tank (7) Pulsation damper 
(4) Pump (8) Centre feed 


marked effect on conductivity. Also at high concentrations 
the conductivity bridge became insensitive because of 
the presence of stray potentials and earthing problems. 
Difficulty was experienced in circulating the small volume 
of solution through the centre compartment, and a special 
roller pump (operating by squeezing the liquid through a 
rubber tube) was developed for the purpose. 

The final experimental method called for very short 
runs — 10-20 minutes. A constant current of four amps 


(3-72 amps correct) was used rather than constant voltage, 
because of its greater significance in determining rate of 
transfer, and because it facilitated interpretation of 
results. The centre compartment flow rate was also kept 
approximately constant at 40 ml/min. (It has been 
found (author, unpublished data), that the rate of material 
transfer through the membrane is independent of flow 
rate in a cell of this type where the Reynolds number 
in the compartments is greater than about 500.) 

(Il) Procedure. The electrode streams were started 
and maintained at 1000 ml/min (500 ml per electrode). 
The centre compartment system was rinsed with distilled 
water, purged with compressed air, and 100 ml of solution 
of approximately correct concentration was added through 
the feed opening. After circulating for about two minutes 
a 25-ml sample was withdrawn for analysis, and the power 
was switched on. At the end of the run, the power was 
cut, and the whole contents of the centre stream were 
blown into a suitable receiver. 

Analysis of the centre stream was achiéved by evaporat- 
ing a 3-ml sample in a deep glass weighing bottle in an 
oven at 150°C, drying for a further four hours and weighing. 

A total of 133 runs was carried out in the cell, covering 
a range of concentration in both streams of about 3-29 
per cent (wt./vol.). The runs were arbitrarily grouped by 
maintaining a fixed electrode stream concentration for a 
series of about twelve centre-compartment strengths. 
1-28 were conductivity monitored and were rejected. 
Runs 29-50 and 78-90 were also rejected because of leakage. 
Both membranes were replaced after run 90, with no 
apparent effect on the results. Because of the nature 
of the equipment, it was inconvenient to maintain a 
controlled temperature, but at all times the operating 
temperatures fell within the range 23-26°C. The short 
running time for each test and the rather complex nature 
of the centre-stream system made impractical direct 
measurement of volume change due to solvent transfer. 

(iii) Results. Cell-current efficiencies, and efficiencies 
calculated from transport-number data are shown in Figs. 
7, 8, 9 and 10. 


3. Discussion 

(a) Transport numbers 

The results show that, for a given arithmetic 
mean concentration, a high concentration ratio 
will result in a correspondingly increased trans- 
port number for the selected ion within the 
membrane. This indicates that, except at low 
concentration differences (say up to 2/1 ratio), 
the concentration gradient between the membrane 
faces is not linear. Replotting Figs. 1 and 2 in 
terms of geometric means or even mean activities 
does not result in a single relationship for the 
experimental data. 
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Since the transport numbers are measured 
indirectly in terms of a balanced potential value, 
they may be considered independent of any 
net electro kinetic transfer across the membrane. 
However with high concentration differences the 
osmotic driving force tending to cause solvent 
transfer will be greater, and, even considering 
the very low hydraulic permeability of these 
commercial membranes, it would be possible for 
solvent ingress (independent of that associated 
with the charged layer near the active resin 
groups) to result in a decreased equilibrium 
concentration within the membrane. As can be 
seen from the figures, a measurement of trans- 
port numbers at a concentration ratio of 2/1, 
as has been accepted practice for most 
investigators, is insignificantly different from the 
ideal condition developed by extrapolation. 

The precautions taken to obtain these highly 
reproducible results have been outlined by 
SoLLNER and Grecor [16]. 


(b) Three-compartment cell 

Efficiencies obtained with the three-compart- 
ment cell were lower than indicated by transport 
numbers at all the concentrations studied. 
This resulted from material transfer by 
mechanisms other than simple electrodialysis. 
Electro-osmosis, which would cause solvent 
transfer from the centre to the electrode com- 
partments, should cause an apparent efficiency 
decrease. Osmosis and diffusion may either 
improve or impair current efficiency, depending 
on the relative concentrations in the centre and 
outer compartments. It is clear, that these 
modes of transfer are not independent of one 
another, and simplifying assumptions could lead 
to wrong conclusions especially when high 
ambient and interstitial concentrations occur. 

However, for each set of runs there is a situation 
where the concentrations on both sides of the 
membrane are equal, and here it may be 
assumed that electro-osmosis is the main factor 
causing discrepancy between cell efficiencies and 
those calculated from membrane transport 
numbers. 

From the difference in efficiencies at these 
points the solvent transfer in millilitres per 


equivalent of NaCl has been calculated, and is 
plotted against ambient concentration in Fig. 11. 
Although the results show considerable scatter 
it is clear that water transfer is greater at high 
concentrations, because of the increased solvent- 
ion interaction. Electro-endosmotic effects are 
probably secondary under these conditions 
because of the high solution conductivity and 
relatively low contribution of the charged 
exchange groups in the membrane pore walls. 
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If the experimental data in Figs. 8, 9 and 10 
are checked using the solvent transfer data of 
Fig. 11 (assuming that for unequal ambient 
concentrations, an arithmetic mean concentration 
is applicable) a reasonably good agreement is 
obtained, showing that in a dynamic system the 
contribution of osmosis and diffusion is small. 


4. CONCLUSION 


It has been shown that the ion-permeable 
membranes “ Permaplex C-10”’ and Permaplex 
A-10”’ will operate selectively in all concentra- 
tions of sodium chloride. For a given arithmetic 
concentration, the membranes show greater 
selectivity as the concentration gradient through 
the membrane is increased. 

A three-compartment cell has been shown to 
give reproducible current efficiency data for the 
membranes over a range of concentrations and 
concentration gradients. 

The results of this work are being confirmed in 
a multi-compartment electrodialysis cell. 
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NOTATION 
a* = Activity of solute. 


E, = Ideal membrane potential, volts. 

E = Measured membrane potential, volts. 

F = The Faraday. 

M = Solution concentration, molar. 
M* = Monovalent cation. 
m* = Monovalent cation passing anion selective 

membrane. 

R = Gas constant, joules deg"! mole~'. 

T = Absolute temperature. 

= Transport number of cation in“ Permaplex C-10, 


1, = Transport number of anion in “ Permaplex A-10.” 
t. = Transport number of anion in “ Permaplex C-10.” 
= Transport number of cation in Permaplex A-10.” 
X~ = Monovalent anion. 


= Monovalent§ anion passing cation selective 
membrane. 
» = Current efficiency, equivalents per Faraday. 
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The kinetic implications of an empirically fitted yield surface for the 
vapour-phase oxidation of naphthalene to phthalic anhydride 


P. H. Prxcuseck* 
The Coal Tar Research Association, Oxford Road, Gomersal, nr. Leeds 


Abstract—A previous investigation resulted in the fitting of an empirical yield surface for 
phthalic anhydride formed in the vapour-phase oxidation of naphthalene. The significance 
of various aspects of this surface are now considered in relation to the reaction kinetics. The 
shape of the fitted surface is shown to exclude all but fractional order kinetics and to provide 
strong evidence for a particular kinetic route. Subsidiary experiments confirm an indication 
that a change in mechanism occurs near to 400°C. The reaction is found to be controlled by 
surface kinetics below this temperature. 


Résumé— Des recherches antérieures ont consisté a ajuster empiriquement la surface de réaction 
VOL. necéssaire 4 la production d’anhydride phtalique par oxydation du naphtaléne en phase vapeur. 
L‘auteur considére ici la signification des différents aspects de cette surface en rapport avec 
la cinétique de la réaction. 
La forme de la surface exclut toute cinétique d'un ordre autre que fractionnaire et semblerait 


indiquer une cinétique particuliére. Des expériences complémentaires confirment une variation 
dans le mecanisme de la réaction autour de 400°C. En-dessous de cette température, la réaction 
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est controlée par la cinétique de surface. 


INTRODUCTION 
In a previous paper [1] a systematic statistical 
study of the vapour phase oxidation of naph- 
thalene to phthalic anhydride was carried out. 
For a particular V,0, K, 50, SiO, catalyst it 
was demonstrated that the phthalic anhydride 
yield surface was well approximated by the 
quadratic expression 
= — 234-89 + 151-61X, + 258-91X, 

+ 810-76X, — 19-95X," 

— 60-86X,,” — 73-07X,” — 50-84X, X, 

— 71-46X, X, — 143-19X, X, (1) 


This expression was found to represent a 
family of hyperboloids as shown in Fig. 1. 
Cross sections of the surface at constant tempera- 
ture are families of ellipses and Fig. 2 shows 
the cross section at 405°C with 70% and 80% 
mole conversion contours. Differentiation of 
equation (1) reveals the existence of a stationary 
point at 401°C in the form of a col or minimax 


(i.e. a point which represents a maximum for 
two variables and a minimum with respect to 
the third). Further calculation shows that the 
maxima at various temperatures increase slowly 
in value both above and below the minimax 
temperature. 

Box [2] has recently drawn attention to the 
fact that “ natural”’ variables, such as tempera- 
ture, pressure and concentration, are used in 
chemical experimentation because of the con- 
venience of measurement. The fundamental 
variables controlling chemical reactions (e.g. 
frequency of molecular collision) are usually 
functions of two or more of the natural variables. 
As a result of this, a yield surface fitted in a 
factor space of natural variables usually contains 
ridge systems. These ridge systems in the 
empirical surface may be related to the basic 
mechanism of the reaction and can be used to 
differentiate between various theoretically 
possible mechanisms [3]. 

Iorre and SHerman [4] have studied the 


*Present address : United Coke and Chemicals Co. Ltd., Sheffield. 
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70% mole conversion The same authors [5] have also investigated 
the steps occurring in the oxidation reaction 


kinetics of the vapour phase oxidation of 
naphthalene and conclude that the reaction rate 
is controlled by diffusion of reactants in the 
catalyst pores at temperatures as low as 380°C. 


and have proposed the following system, 


4 NAPHTHAQUINONE 


ks cCH— CO 


In the course of an investigation into catalysts 
for the oxidation of SO, to SO,, Tanpy [6] 
reported that an approximate melting point of 
400°C for mixtures of V,0O, and K,SO,. The 
catalyst used in obtaining equation (1) contained 
V,0, and K,SO, in similar proportions to that 
used by Tanpy. 


Tue Fittrev SURFACE 


Consideration of the form of the fitted yield 
surface (equation 1) reveals several significant 
facts ; 

1. The presence of a ‘ minimax’ within the 
experimental region raises the possibility 
of a change of kinetic mechanism around 
this temperature. 


The attenuated ridge systems occurring 
indicate a considerable degree of factor 
inter-dependence. 

Although the empirical surface was fitted 
to experimental results spread over 
unusually wide regions of factor space and 
included mole conversions between 60 and 
85% the deviations of the experimental 
results from the predicted results are 
little greater than those expected from 
the known experimental errors. The true 
yield surface is, therefore, very similar 
in form to the model obtained. 
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The kinetic implications of an empirically fitted yield surface for the vapour-phase 


Table 1. 


40/80 Mesh 
Property sample sample 


100/150 Mesh 


Average particle diameter cm 


External surface area of one average 
particle em? 


Average pore length cm 
Total surface area m?/g B.E.T. method 38-9 


Percentage pore volume 


Particles per gram 


External surface area m?/g 


26 x 


21-1 x 


0-6 x 10° 


1-25 x 10°? 


491 x 10% 


03 x 10% 


It seemed desirable at this stage to investigate 
two main aspects of the yield surface : first to 
determine whether the reaction taking place 
was controlled by diffusion in the catalyst pores 
and whether any change in mechanism occurred 
near to the minimax temperature ; secondly, 
by adopting the kinetic system proposed by 
lorre and SuerMan, to compare the fitted surface 
with surfaces generated by various theoretically 
derived rate equations, 


REACTION RaTE AND PARTICLE 


If the reaction rate is retarded by diffusion of 
reactants in the catalyst pores then an increase 
in the number of pores available for reaction 
should give a corresponding increase in reaction 
rate. Accordingly a sample of microspherical 
V,0,/K,SO, on silica gel catalyst, similar to 
that used in the previous work, was prepared. 
This catalyst was accurately sieved between 
40 and 80 mesh screens and divided into two 
parts. One part was then ground down and a 
sample sieved out between 100 and 150 mesh. 
Various physical properties of these two samples 
were determined and are recorded in Table 1. 

As the two samples are from the same original 
batch the number of pore mouths per unit surface 
area will be almost identical. Thus the number 
of pore mouths available in a given weight of 


catalyst sample is proportional to the external 
surface area and the finer material should contain 
1.73 times as many as the coarser sample. As, 
however, the 100-150 mesh particles are not 
spherical the external surface area will be greater 
than that calculated and the factor is likely 
to be somewhat greater than 1-73. 

Using identical weights of these samples, the 
rate of oxidation of naphthalene was examined 
under the same operating conditions, chosen 
in such a way as to obtain incomplete conversion 
of the naphthalene. The experimental apparatus 
and analytical techniques used in this work 
were the same as those described previously [1]. 
Fig. 3 shows the Arrhenius plot of first order 
rate constant with temperature. The energies of 
activation obtained from the slopes of the 
straight line portions were, 


40/80 mesh catalyst 26-0 kcals/mole 
100/150 mesh catalyst 28-0 kceals/mole 


These may be compared with the results 
obtained by CaLDERBANK [7] (26 kcals/mole), 
Mars and vAN Krevecen [8] (22 kcals/mole) 
and Iorre and SHerMan [4] (27-4 kcals/mole). 

It is clear that the expected increase in rate 
has not taken place as the constant obtained 
below 400°C for the finer catalyst does not 
exceed that from the coarser material by more 
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Remarks 
Calculated assuming 
perfectly spherical 
particle 
38-9 Not determined 
separately 
34 35 
4-42 x 104 82:7 x 104 
| 0-0093 0-0161 Ratio 1-73 


— 
Reciprocal absolute temperoture 10° 
Fic, 3. 


than a factor of 1-15. Also the energy of activa- 
tion obtained shows that the Arrhenius relation 
is obeyed up to 400°C and consequently, as 
shown by WHee er [4a], surface kinetics must 
be rate controlling up to this point. Between 
400° and 500°C the difference between the two 
catalysts is more marked and it is possible that 
diffusion plays a greater part in this region. 
The results in the region of 500°C were the 
least reliable, and this point would require more 
careful investigation. There would seem to be 
no doubt, however, that the kinetic mechanism 
undergoes modification around 400°C as indicated 
by the rapid change of slope in Fig. 3. This 
provides support for the contention that the 
occurrence of a *‘ minimax’ in the fitted surface 
is not a chance effect but is the result of fitting a 
quadratic expression to a true surface containing 
a major axis which has a change of slope at this 
point. 
Kinetic Orper THE YIELD SURFACE 
Using as a model, the kinetic system proposed 
by Iorre and SuHerMan [5] various cases may be 
considered and integrated rate equations deduced. 


Case 1 
Disappearance of naphthalene half order, 
subsequent steps first order. 


k t 
Case 2 
All steps in reaction first order. 
ky ks 
(kg + ky) — (hy + ky) + ke) 


Cy (eo 


e 


+ 
ks — (kg + ky) 


Case 3 


Disappearance of naphthalene second order, 
subsequent steps first order. 


ky ks 


em 
CN, ks (ks + ky) 


ky ky 1 
kes (ky + + (hy + hy) t 


C, 

Cy 
conversion, these rate equations may be made 
equivalent to the fitted equation. The rate 
equations and the fitted equation then contain 
the same variables. Temperature appears in the 
rate equations by the dependence of the rate 
constants, initial naphthalene concentration, 
Cy, and air/naphthalene ratio are inversely 
proportional and contact time appears directly 
in all cases. 

The rate equations may thus be compared 
with empirical equation on the basis of the yield 
surface generated in the same factor space. 
For convenience a two dimensional factor space 
at constant temperature was first considered. 
Using the rate constants given by lorre and 
SHERMAN, values were obtained by solution of 
the rate equations. The surfaces were then 
plotted on the same scales as those used in fitting 
the empirical surface. The sections so obtained 
are shown in Fig. 4. 

It is seen that as the order of the first step 
increases, so the slope of the major axis (in the 
factor space considered) changes from negative 
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to positive. For this system, therefore, orders negative slopes and orders greater than unity 
less than unity will give rise to surfaces having will give positive slopes. 
It is evident that the fractional order equation 
generates a yield surface of strikingly similar 
Case 1 form to that obtained empirically and _ this 
Half order equation surface was considered in greater detail. The 
values of the various rate constants given by 
IorreE and SHERMAN do not predict maximum 
conversion at the same position in the factor 
space as the fitted equation, neither do they 
predict the same level of maximum. This is 
not surprising since the rate constants will vary 
in magnitude and ratio with the composition and 
physical form of the catalyst used. The values 
of the constants at any temperature will 
determine 


(i) The maximum conversion 


(ii) The position(s) in the factor space at 
which maximum conversion is obtained 
(iii) The “ spread” of the contours. 


It was also found that variation of the ratios 
ky ky, and kg/k, caused variations in the slope 
of the major axis of the constant temperature 
section. Consequently, this slope is likely to 
vary with catalyst form and composition. 

Considering the complete three dimensional 
factor space (i.e. including temperature) the slope 
of the entire surface with respect to the tempera- 
ture axis will be controlled only by the rate 
of change of the values of the various rate 
constants with temperature (i.e. dependent upon 
the activation energy of the various steps). 
Within the region where surface kinetics are rate 
controlling this rate of change will be independent 
of catalyst and thus all catalysts should show 
a parallel major axis. In order to check this 
point the position of maximum conversion was 
determined from the fractional order equation 
at various temperatures and Figure 5 shows a 
comparison of the slopes obtained with those 
derived from the empirical equation. Once again 
a large measure of agreement exists. 


6 
3 


CONCLUSIONS 


The comparison of reaction rates obtained over 
different size ranges of the same catalyst indicates 
that, contrary to the results obtained by lorrr. 
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VOL. 
Case 2 
| 20 
Case 3 
Second order equation 
max ~ 
Fic, 4. 


P. H. Pincnspeck 


Major axis slope X,/,, Major axis slope X3/y. 


Empirical 
equation 


Air naphthalene rotic 
(ogerithmic scale) 


Fie. 5. 


and SHERMAN [4], the reaction rate is controlled 
by surface kinetics up to temperature of about 
400°C. In this respect the results support the 
views of Mars and vaN KREVELEN [8] who also 
employed fluidized catalysts of a similar composi- 
tion. From the data of Ilorre and SHerman [4] 
it can be deduced that their catalyst was of 
unusually high apparent density (1-91) and con- 
sequently had a low porosity. It is also well 
known that variations in the treatment of silica 
gel can lead to widely differing pore structures. 
The results of lorre and SHERMAN may, therefore, 
be explained in terms of the physical form of the 
catalyst used. 

The temperature at which the ‘ minimax ’ 
occurs in the fitted surface is shown to be 
associated with a change in kinetic mechanism, 
a change in state of the active catalyst com- 
ponents also occurs near this temperature. 

From the comparison of the types of yield 
surface generated by various rate equations strong 
evidence has been obtained indicating that the 
vapour phase oxidation of naphthalene, under 
high conversion conditions, occurs with fractional 
order kinetics. The direction of the axes of the 
empirical surface would suggest an order between 
0-8 and 0-9. The agreement between the empirical 


and theoretical surfaces provides strong support 


for the kinetic system proposed by lorre and 
SHERMAN [5], and the slope of the surface with 
respect to the temperature axis indicates that the 
values of activation energies deduced by these 
authors for the various steps are reasonable. 
This evidence is more powerful than the normal 
kinetic data in that it represents the results of 
some seventy experiments carried out over wide 
regions of factor space considered as a group. 

The investigation demonstrates that the purely 
empirical approach, when treated systematically, 
can throw considerable light on the basic 
mechanism of the reaction under investigation 
and provides a new approach to complex 
multivariable reactions. 
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NOTATION 


Y = % Mole conversion of naphthalene to phthalic 
anhydride 


X, = log, Air/Naphthalene Ratio litres /y 
X, = log yo (Contact Time x 10) sec 


Xs, = (Temperature — 330) 


Cy, = Initial concentration of naphthalene moles /litre 


C, = Final concentration of phthalic anthydride 
moles /litre 


k, = Specific Rate Constant Naphthalene - 1:4 
Naphthoquinone 

kg = Specific Rate Constant Naphthalene + Phthalic 
anhydride 


kg = Specific Rate Constant 1 : 4 Naphthoquinone > 
Phthalic anhydride 


kg = Specific Rate Constant 1:4 Naphthoquinone 
— Complete combustion 


ks = Specific Rate Constant Phthalic anhydride -> 
Maleic Anhydride 
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Fugacities in simple gas mixtures 
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Division of Chemical Engineering, University of California, Berkeley, California 


( Received 6 November 1956) 


Abstract—A method is presented for rapidly calculating the fugacity of a component in gaseous 
solution using only critical data for the pure components in the solution. The method is based 


on an extension of the theorem of corresponding states proposed by GuGGeNHEIM and McGLasHaNn 
and is applicable to multicomponent mixtures of simple gases at moderate densities. 


Résumé—L auteur présente une méthode de calcul rapide de la fugacité d'un gaz dans un 
mélange de plusieurs gaz 4 partir des données critiques des constituants purs. 
Cette méthode basée sur une généralisation du principe des états correspondants de GUGGEN- 


Fucacities of components in a gaseous mixture 
must frequently be estimated in problems con- 
cerned with physical or chemical equilibria at 
advanced pressures. Since the extensive experi- 
mental data required for a rigorous solution are 
very seldom available it is customary to utilize 
approximating methods; the simplest of these 
is the Lewis fugacity rule which states that 


fi (1) 
where f, = fugacity of component i in the mix- 
ture. 


Sf? = fugacity of pure component i at the 
temperature and pressure of the mix- 
ture. 


y, = mole fraction of component 7 in the 
mixture. 


The Lewis fugacity rule assumes that Amagat’s 
law of additive volumes is valid for the pressure 
range from zero pressure to the pressure of the 
mixture ; this assumption is known to be invalid 
for many if not indeed for most cases. 

More reliable methods for estimating fugacities 
in mixtures are the method of Jorre [5] who 
utilizes Kay's pseudocritical concept, and the 
method of Repiicu and Kwone [7] who use for 
the mixture a two-parameter equation of state 
whose constants are a simple algebraic combina- 
tion of the constants for the pure components. 


HEIM et McGLASHAN peut étre appliquée a des mélanges de gaz simples a densités moyennes. 


The method proposed in this paper does not 
assume Kay’s rule nor the validity of any parti- 
cular equation of state. Instead, it is based on 
an extension of the theorem of corresponding 
states proposed by GUGGENHEIM and McGLasHaN 
[2]. This method is very easy to use and at 
moderate densities may be expected to give good 
results for simple mixtures ; that is, for mixtures 
of nonpolar molecules whose force fields are 
spherically symmetric. 

For a simple binary mixture the second virial 
coefficient B,,, is given by a quadratic function 
of the composition (6) : 


= + 2 By + By,’ (2) 


where B, and B, are the second virial coefficients 
of the pure components 1 and 2 respectively, and 
where B, , is a constant for the binary mixture. 
and McGuiasnan [2] have shown 
that the theorem of corresponding states can be 
used not only to predict B, and B, but also B, , 
provided the critical volume and critical tem- 
perature of both components are known. 

The second virial coefficients for the pure 
components are given by 
B, T 

where T = temperature ; V = volume per mole ; 
and subscript ¢ denotes critical. 


(3) 


VOL. 
6 


956/57 


Fugacities in simple gas mixtures 


The data for numerous simple gases show that 
the function @ is a universal function which is 
very nearly the same for all simple gases. 

GuGGENHEIM and McG.LasHan have examined 
the data for various simple binary mixtures and 
have shown that the coefficient B,, can be 


(4) 


where the function @ is the same as that for the 
pure components. The terms V.,, and 7T,, , are 
found from the corresponding terms for the pure 
components by means of the common mixing 
rules as recommended by Hirscure.per et al., 
in their work on transport properties [4] : 


Vey, = 4 + (5) 


In most problems it is more convenient to use 
the compressibility factor rather than virial 
coefficients. The relationship between the com- 
pressibility factor z and the virial coefficient B 


If higher terms in the virial expansion are 
neglected, then at constant pressure and tem- 
perature z is a linear function of B and from 
equation (2) the relationship between z,,,, the 
compressibility factor for the mixture, and the 
composition is given by 


= + 2 + (8) 


where z, and z, are the compressibility factors for 
pure components 1 and 2 respectively, at the 
temperature and pressure of the mixture. (If it 
is assumed that z, , = $(z, + 2,) then equation 
(8) reduces to Amagat’s law of additive volumes.) 

The compressibility factors for the pure com- 
ponents are readily found from the generalized 
compressibility charts if the critical pressures and 
temperaturse of the pure components are avail- 
able. To find the compressibility factor z, , it is 
necessary to define a critical pressure P., . which 
can be found by extending the previously men- 
tioned mixing rules as follows : 


The critical pressure P, for a pure component 
is given by 
2, RT, 
(9) 
where z, = compressibility factor at the critical 
point. For P., . we therefore write 


(10) 


where V., , and 7., , are given by equations (5) 
and (6). The term z, , is some average between 
Ze, and %,; fortunately z, does not vary much 
for simple molecules and therefore any method 
for averaging will give very nearly the same 
result as any other method. It is simplest to 
adopt the arithmetic average and equation (10) 
then becomes 


+ 
The fugacity of a component in a binary mixture 
can readily be calculated from equation (8) and 
can be expressed in the form 


f= nt (12) 

where y, = activity coefficient of component 1 in 
the mixture. 

f,° = fugacity of pure component 1 at the 


same temperature and pressure as 
those of the mixture. 


The ratio of the fugacity of a component in a 
binary mixture to its partial pressure is given by 
an equation previously derived by Jorre [5] : 


The ratio of the fugacity of a pure gas to the 
total pressure is related to the compressibility 
factor of that gas by 


1—% 
In(f°/P), = —| 4P (14) 
and similarly for the mixture as a whole 


In(f/P) nix = (15) 
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Substituting equation (15) into equation (13), 
inverting the order of integration and differentia- 
tion, and then substituting equations (8), (14), 
and (12) yields a simple result for the activity 
coefficient y : 

In y, = Ky,” (16) 
where 


K = —In(f°/P), (17) 
Similarly 
and 


So = (18) 
In y, = Ky,* (19) 


The constant K in equation (19) is the same as 
that in equation (16). 

The three terms which comprise K can readily 
be found from the generalized fugacity coefficient 
charts which express f/P as a function of reduced 
pressure and reduced temperature. The critical 
constants to be used in the evaluation of the first 
term are given by equations (11) and (6). 

If Amagat’s law holds for the entire pressure 
range from zero pressure up to the pressure of 
interest, then K = 0, and y, = y, = 1. 


APPLICABILITY 


In the derivation given above all virial co- 
efficients higher than the second were neglected ; 
the proposed method of calculation therefore is 
only valid for moderate densities. Equation (8) 
predicts that the volume change in mixing is a 
symmetric parabolic function of the composition 
which is in agreement with experimental results 
for simple mixtures at moderate densities. The 
data for the nitrogen-ethylene system at 50°C (8), 
for example, show that this symmetric parabolic 
relationship holds quite well for pressures up to 
about 100 atmospheres. Accoridng to Amagat’s 
law, of course, the volume change on mixing is 
zero. 

Some results for the nitrogen-ethylene system 
are shown in Table 1. The activity coefficients 
predicted from equation (16) and the generalized 
charts are compared with those calculated by 
Bennett [1] from the data of Hacensacn and 


Comincs [8]. The agreement is good. 


Table 1. Activity Coefficients for Nitrogen- Ethylene 
Miztures at 50°C 


YN, 


Cale. | Ref.{1]| Cale. | Ref. [1] 


The approximate range of temperature and 
pressure for which this method is applicable is 
that for which virial coefficients higher than the 
second can be neglected without serious error ; 
this range is for densities up to about one-half 
of the critical. As a rough rule, therefore, the 
method is applicable for the range where 


Vv 


Substituting V = zRT,/P and using a value of 
0-29 for the compressibility factor at the critical 
point the valid range is approximately given by 


p 

—— > 0-58 

Pr ? 

where the critical properties of the mixture may 
be estimated by Kay’s rule. 


EXTENSION TO MULTICOMPONENT 
MIXTURES 
The above derived method for a binary mixture 
can readily be extended to a multicomponent 
system by the well-known methods summarized 
by Wout [8]. Since equations (16) and (19) for 
the binary case are the result of a two-suffix 
expansion for the excess free energy, extension 


1l4 


Vol 
6 


195¢ 


| 
Mol 
Pressure| Fraction — 
(atm) | Nitrogen 
50 0-0 1-18 1-19 1-00 1-00 
0-2 1-10 1-12 1-01 1-01 
0-4 1-06 1-06 1-02 1-02 
0-6 1-02 1-02 1-06 1-07 
0-8 1-01 1-01 1-10 
10 1-00 1-00 1-18 1-18 
100 0-0 1-41 1-45 1-00 1-00 
0-2 1-25 1-27 1-01 1-01 
lll 1-14 1-06 1-06 
0-6 1-06 1-06 Ill 1-14 
0-8 1-01 1-01 1-25 1-27 a 
10 1-00 1-00 141 145 
|| 
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to a multicomponent system requires only con- NOMENCLATURE 
stants for all the binaries in the multicomponent —_g ~ second virial coefficient 
mixture. The activity coefficient for a species r 7, — fugacity of component i in the mixture 


in a mixture of n components if given by J®° = fugacity of pure component i at the temperature and 
pressure of the mixture 


n n n 
K = defined by equation (17) 
= >) D Ky (20) P = pressure 
i=1 i=l j=l R = gas constant 


The coefficients Ky are evaluated from the T = temperature 
generalized fugacity chart as indicated previously; 


. y = mol fraction 

z = compressibility factor 
K,;, = 2ln(f/P),; — In P),; — In P y = activity coefficient 
ij (f/P)s (f°/P); (f°/P); 

When i = j the coefficient K vanishes as expected. Subscript R denotes reduced 
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Prediction of ternary liquid-liquid equilibrium from binary data 
(for the system furfural—iso-octane—benzene) 
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Abstract—The distribution of benzene in the ternary system furfural — iso-octane — benzene 
was predicted from binary data, obtained by total pressure measurements on the systems 
benzene - iso-octane, benzene -furfural and iso-octane-furfural. The distribution was 
predicted using the two - suffix equations of van Laar [2] and the semiempirical methods of 
Hicpesranp [5] and and Frrepianp 

The system carbon tetrachloride — acetic acid was also investigated by the total pressure 
method. 

The binary data obtained from these total vapour pressure measurements can be taken VOL 
as fairly reliable, since azeotropes were picked up quite accurately and Repuicu and Kiser plots 6 
showed that the results were thermodynamically consistent. The method of calculation ensures 

that the final activity coefficients are a- very good representation of the original pressure 

measurements. 

This binary data, when applied in the van Laar equations for the prediction of ternary 
equilibrium, gives a distribution which is consistently 1-2 mole % of solute low in the iso-octane- 
rich phase when compared with Henry's experimental distribution. The distribution predicted 
by Hi-peBranp’s method is more accurate than the van Laar prediction in the region near 
the plait point, but is fallacious in the near binary regions in that it indicates a molar solutrope. 

The ternary van Laar distribution obtained using the binary data and the binodal curve 
is accurate in the region near the plait point, but remains 1-2 mole % of solute in error near 
the binaries. 

Scneret and Frrepianp’s correlation is generally more satisfactory than van Laar’s 
except near the binaries, where it indicates a molar solutrope, which does not in fact exist. 


Résumé— L'auteur étudie la distribution du benzéne dans le systéme ternaire furfural-isooctane— 
benzéne en fonction des pressions totales des systémes binaires : benzéne—isooctane, benzéne— 
furfural, isooctane—furfural. I) étudie cette distribution en utilisant les équations de Van Laar 
et les méthodes semi-empiriques de HiLDEBRAND et de SCHEIBEL et FRIEDLAND. 

Le systéme tétrachlorure de carbone—acid acétique a été également étudié parla mesure des 
pressions totales. 

On peut admettre que les données binaires obtenues a partir des mesures de pression sont a 
peu pres exactes, étant donné que les azéotropes sont retrouvés avec exactitude et que les courbes 
de Repuicu et Kister donnent des résultats compatibles thermodynamiquement. La méthode 
de calcul garantit que les coefficients d'activité sont en accord avec les pressions initiales. 

L’auteur applique ces données binaires a l'équation de Van Laar en vue d’obtenir l’équilibre 
ternaire : il en résulte une distribution inférieure de 1 4 2 mole % de soluté dans la phase riche 
en isooctane comparée a la distribution expérimentale de Henry. La distribution prévue par la 
méthode de HitpeBRAND est plus précise que celle de Van Laar dans la région du point de 
courbure, mais ne convient pas dans les régions binaires en ce sens qu'elle indiquerait un solutrope 
molaire. 

La distribution de Van Laan basée sur cette méthode, données binaires et courbes binodales, 
est satisfaisante dans le voisinage du point de courbure, mais accuse une erreur de 1 & 2 mole %, 
de soluté prés des zones binaires. 

La relation de ScuriBe. et FraiepLanp est en général plus satisfaisante que celle de vaN Laar 
sauf dans les régions binaires ow elle semblerait indiquer un solutrope molaire qui, en fait, n'existe 
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INTRODUCTION 


Tue aim of the work was to investigate the 
accuracy with which certain ternary liquid- 
liquid equilibrium relationships could be predicted 
from binary data and some ternary data. Various 
mathematical procedures have been proposed 
by Wout [1], van Laar [2], Benepicr [3], 
Rep.icu and Kister [4] and others. In addition 
to these, empirical or semiempirical methods 
have been suggested by HitpeBranp [5] and 
ScHEIBEL and FriepLanp [6]. These methods 
must be reliable since it is essential that the 
distribution of one liquid between the other two 
solvents be known accurately before the process 
can be applied to industry. 

Binary data were obtained by the measurement 
of total vapour pressures on a range of binary 
liquid mixtures. From these readings the activity 
coefficient-composition curves were calculated 
by the method due to Barker [7]. 

The ternary distribution was predicted by 
combining the data for the three relevant binaries. 
Experimental ternary data were supplied by 
C. J. Henry who is studying at Sydney a quater- 
nary system with iso-octane — furfural — benzene 
as one ternary. 


APPARATUS AND MATERIALS 


A sketch of the total pressure apparatus is 
shown in Fig. 1. The vapour space was 
warmed by an isotape to prevent distillation and 
consequent alteration in liquid mixture composi- 
tion. Provision was made for agitation under 
vacuum by means of a nickel paddle actuated 
by the rotating magnet of a magnetic stirrer. 

A Cenco-Megavac pump was used to obtain 
a vacuum of the order of 0-1 microns of mercury. 
Water at 25-0°C for the bath in which the phial 
was immersed was circulated from a large 
laboratory tank by means of a centrifugal pump. 
Pressure differentials could be read to 0-01 mm 
of mercury with a cathetometer. Liquid nitrogen 
in vacuum flasks was used in the vapour traps 
and for freezing purposes. 

Components of a very high degree of purity 
were supplied by R. G. H. Prince and 
C. J. Henry, who were engaged on liquid-liquid 
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Fic. 1. Sketch of Total Pressure Apparatus. 


equilibrium researches in the Chemical Engineer- 
ing Department of the University of Sydney. 
Important physical constants for these com- 
ponents are listed in Table 1. 


Table 1 
Carbon | Acetic 

Properties | Benzene | Isooctane Furf tet. acid 
Boiling 80-1°C 161-7°C | 76-7°C: | 117-9°C 
point 
Density, 08736 | 0-68789 1-15482 | 15846 11-0438 
25°C g/ml ig/ml jg/ml g/ml 
Refractive 1-38896 
index, np 
Freezing | 5-47°C 38-8°C 16-55°C 
point 

These components were stored in sealed 


Winchesters and were carefully dispensed from 
automatic burettes fitted with drying tubes. 
In this way the liquids retained their high degree 
of purity throughout the experiments. 


mm bore vacuum 
capillary pump 
x 
a7 | all 
quickfit 
joint | 
| 824 quickfit tube 
Air trap 
VOL. | 
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EXPERIMENTAL PROCEDURE 

The binary mixtures were calculated roughly 
in terms of the volumes of the pure components. 
These volumes were introduced into a stoppered 
50 ml conical flask which was weighed accurately 
after each addition of liquid. 

The vapour pressure was determined in the 
following fashion:— the phial was almost 
filled with the liquid mixture, which was then 
stirred gently to remove air bubbles and frozen. 
The apparatus was evacuated, the mercury 
raised into the U-tube and zero levels read. 
The mixture was then thawed, stirred, and 
allowed to come to equilibrium at 25-0°C, when 
final levels were read. The change in mercury 
levels in the U-tube was the apparent vapour 
pressure of the mixture in the phial. Actually 
this figure ‘had to be corrected for the presence 
of occluded air and other permanent gases. 
These can only be removed by vacuum sublima- 
tion and this alters the liquid mixture composi- 
tion, a feature which was undesirable in this 
work. 


Table 2. System Carbon Tetrachloride — Acetic Acid 


Run no.| Mole fraction | Mole fraction | Total vapour 
carbon 
tetrachloride i i mm Hg 


the vapour volume slightly and noting the change 
in total pressure. Knowing the initial and final 
volumes from the dimensions of the apparatus, 
the pressure of permanent gases causing such 
pressure change was then calculated, from 
Boyle’s Law. This pressure for one particular 
vapour volume was then subtracted from the 
total pressure at that volume, giving the mixture 
vapour pressure. 


EXPERIMENTAL READINGS 


The literature [8] gives the minimum boiling 
azcotrope at 0-927 mole fraction carbon 
tetrachloride. 


Table 3. System Benzene — Iso-octane 


Total vapour 
pressure 
mm Hg 


Run no. | Mole fraction | Mole fraction 
benzene 


0-0000 


1 
2 
3 
4 
5 
6 
7 
8 
9 


0-0000 
0-0209 
0-0386 
0-2987 
0-5030 


The apparent vapour pressures of the pure 
components were in error by from about 1 mm 


of mercury for those with low vapour pressures 
to about 4mm of mercury for those with high 
vapour pressures. Intermediate mixtures were 
corrected proportionately to their compositions. 
These corrections were also checked by altering 


The literature [8] gives the azeotrope as 
0-981 mole fraction benzene. 


Table 4. System Benzene — Furfural 


Run no.| Mole fraction | Mole fraction | Total vapour 
furfural pressure 
tp mm Hg 


1-0000 146 
0-9479 
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| 
1-0000 49-0 
0-0520 0-9480 59-8 
0-1036 0-8964 66-0 
0-2980 0-7020 80-9 
0-6930 0-3070 94-8 
0-9000 0; 1000 98-7 
0-9468 0-0532 99-2 
1.0000 0-0000 95-2 
1-0000 15-2 
09614 83-0 
0-7018 87-0 
0-4970 99-2 
0-7000 0-3000 110-0 
0-9180 0-0820 121-0 
0-9640 0-0360 124-0 
0-9780 0-0220 118-6 
1 1-0000 0-0000 1110 
0-0000 
0-0521 
0-1027 0-8973 21-0 
0-3045 0-6955 47-4 
0-5060 0-4940 64-6 
0-7035 0-2965 78-0 
0-8997 0-1003 90-2 
0-9526 0-0474 92-7 
1-0000 0-0000 95-2 


Table 5. System Iso-octane — Furfural 
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0-2702¢ + 1-588 272 + 1-096 


The terminal solubilities as determined by Mr. Henry 
are : 

0-0431 mole fraction iso-octane 
(1) 

0-9569 mole fraction furfural 

0-9461 mole fraction iso-octane 


(2) 
0-0539 mole fraction furfural 


Total pressure measurements were also made 
on six ternary equilibrium mixtures. The 
experimental readings are presented in Table 6, 
where they are compared with the pressures 
calculated from the formula 


P= + Po + Ps 


where the activity coefficients are calculated 
by the van Laar equations. 


EXPERIMENTAL RESULTS 

Complete activity coefficient-composition rela- 
tionships for three of the four binaries studied 
are given by Figs. 2-4, inclusive. The extrapo- 
lated end values, which are the constants used 
in the vAN Laar equations etc., are marked. 
It should be noted that these have to be divided 
by 2-303 to reduce them to the common logarithm 
basis on which the equations are worked. 

For the almost immiscible system iso-octane— 
furfural these end values were calculated from 
mutual solubility data. 

The two-suffix van Laar equations for the 
activity coefficients in the ternary reduced to : 


1-308 27+ 0-5842¢ + 1-552 ap 2% 
+ 0-986 ry +,1-041 ac)? 


log yz = 
0-678 1-068 x 
Run no.| Mole fraction | Mole fraction | Total vapour wohegitieaid od 
iso-octane Surfural pressure log yc = 0-474a 2 + 0-871 27 — 004027, 
a4 on + 0-961 + 1-476 25) 
1 0-0000 1-0000 1-6 
2 0-0170 0-9830 26-0 
8 0-0450 0-9550 52-2 28 o 
4 0-9550 0-0450 54-0 
5 09790 0-0210 52-0 
6 0-0000 49-0 


wn 
O-6F 


02 04 


Fic. 2. Plot In y vs. 2 for carbon tetrachloride — acetic acid. 
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Mole fraction benzene - xc 


Fic. 3. Plot In y vs. 2 for benzene — iso-octane. 
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Mole fraction benzene - xc 


Fie. 4. Plot In y vs. 2 for benzene — furfural. 


Figures 5-8, inclusive, present a comparison 
of ternary equilibrium data predicted by various 
mathematical and empirical methods with 
Henty’s experimental data. 


Discussion OF EXPERIMENTAL ResvuLts 
(A) Binary Data 

The theory underlying the method of calcula- 
tion of the binary activity coefficient-composition 
curves is such that the final results are a very 
good representation of the original total pressure 
measurements. The constants obtained by 
extrapolating the curves were checked by 
applying the van Laar equation to the experi- 
mental points nearest 2, = 0 and 2, = 1-0 and 
the agreement between extrapolated and 
calculated constants was always within 1% of 
error. 

These constants should be quite reliable, 
since such azeotropes as existed were picked up 
quite accurately and Repiicu and Kiser plots 
for the three miscible binaries gave nett areas 
under the curves = zero with an error of less 
than 1%, proving the data to be thermodyna- 
mically consistent. 


(B) Prediction of Ternary Data 
(a) From binary data only 
Figs. 5 and 6 show the degree of agreement 


between experiment and prediction using binary 
data alone. From Fig. 5 it can be seen that 


Benzene (C) 


Fic. 5. Prediction vs. Experiment for system iso-octane - 
furfural —-benzene (triangular). 
Key 
Predicted van Laan results — upper curve. 
———- Henry’s experimental data. 


——Henty Z 
Van Loar 
— — Hildebrand 
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xc 
Fic. 6. Prediction vs. Experiment for molar distribution 


in system iso-octane -furfural—benzene from binary 
data only. 
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Prediction of ternary liquid-liquid equilibrium from binary data 


the two-suffix van Laar equations predict the 
binodal curve quite accurately in the region near 
the binaries, but fall down near the plait point. 
The prediction of the tie lines is far from accurate 
and Fig. 6 shows that the predicted distribu- 
tion is consistently 1-2 mole % of benzene too low 
in the iso-octane rich phase. 

The tie lines resulting from van Laar predic- 
tions are almost horizontal in the lower region of 
the triangular diagram and as the percentage of 
solute increases the slope of the tie lines towards 
the furfural apex increases. At no stage is the 
slope as great as that of Henry’s experimental 
tie lines and the discrepancy between tie line 
extremities is sometimes as much as 4 mole 
% of benzene. 

The van Laar distribution giving almost 
horizontal tie lines comes very close to predicting 
a molar solutrope, while the H1LDEBRAND corre- 
lation does in fact predict one at 13 mole °%, ben- 
zene in each phase. No such solutrope exists in the 
system and although the HitpEBRAND method 
correlates well near the plait point, it can in 
general be said that these binary predictions are 
not sufficiently acucrate for application to the 
design of industrial extraction equipment. 


(b) From binary data and some ternary data 

Figs. 7 and 8 show the improvement in the 
accuracy of prediction gained by using Henry’s 
experimental binodal curve in conjunction with 
the binary constants obtained from total pressure 
measurements. 

The van Laan prediction is still consistently 
1-2 mole % of solute in error in the lower portions 
of the diagram, but near the plait point, i.e. 
at about 30 mole % benzene, this error falls to 
0-5 mole % of solute. From a theoretical view- 
point this improvement is pleasing, but from a 
practical point of view it is of little importance, 
since it is generally the lower section of the 
ternary diagram which is used in the design of 
extraction processes. 

The and Friep.Lanp correlation, like 
that of HILDEBRAND, predicts a molar solutrope at 
about 10 mole % benzene and is in error by 
as much as 3 mole % of benzene in either phase 
at concentrations lower than 10 mole % benzene. 


Benzene (C) 


Fic. 7. Equilibrium tie lines for the system iso-octan: — 
furfural — benzene (triangular). 
Key 


Predicted from binary data given ternary satura- 
tion curve. 
— Henry's experimental tie lines 
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Fic. 8. Prediction v. Experiment for molar distribution 
is system iso-octane — furfural — benzene — given binodal 
curve. 


Above this, however, the ScHEerBEL and Friep- 
LAND prediction is by far the most satisfactory 
and is consistently less than 1 mole °%, benzene in 
error. Once again, however, it must be admitted 
that this prediction is unsatisfactory for industrial 
application because of the high errors in the 
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near-binary regions of the ternary distribution 
diagram. 

The improvement in ternary prediction result- 
ing from the use of the three-suffix MarcuLes 
equations (which contain a ternary constant in 
addition to the original binary constants) were 
cursorily investigated. As far as could be judged 
from the few calculations made, the effect of 
introducing the ternary constant was to bring 
the predicted distribution into very much closer 
agreement with the experimental results. The 
most convenient method of estimating the 
ternary constant appears to be from a knowledge 
of one equilibrium tie line. The constant C 
can be solved for directly by equating the 
activities of the solute (expressed in terms 
of mole fractions, binary constants and the 
ternary constant C) at the two extremities of 
the tie line. 

The calculations entailed in fully applying the 
three-suffix ternary MARGULES equations are very 
tedious but it is felt that the increased accuracy 
of prediction would justify the extra work. 


(C) Measurements on Ternary Equilibrium 


Miz ures 

Table 6 
Mole Mole Mole Pressure Pressure 
fraction | fraction | fraction as as 


[2] van Laar J. J. Z. Phys. Chem. 1910 72 723. 


J. W. Kenyw : Prediction of ternary liquid-liquid equilibrium from binary data 


REFERENCES 
Wout Kurt Trans. Amer. Inst. Chem. Engrs. 1946 42 215. 


Table 6 illustrates the excellent agreement 
between experimental and _ predicted total 
pressures for ternary equilibrium mixtures. The 
error in most cases is less than 1%. 

Predicted figures came from applyidg the 
formula 


P = + Po + Ps 


where the activity coefficients are calculated by 
the van Laar equations. 
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NOMENCLATURE 


A = constant in MarGuLes, VAN Laar, equations 
= limit of log y as z > 0 
= component of a solution 

a = activity 

B = component of a solution 

C = constant in 3-suffix ternary MarGuLes equations 
= component of a solution 

P = total vapour pressure, mm mercury 

P, = vapour pressure of component 1, mm mercury 


AB = component A in a B-rich solution, etc. 


[3] Benepict M., Jomnson C. A., Socomon E., and Rupen L, C. Trans. Amer. Inst. Chem. Engrs. 1945 41 371. 
(4] Repurcn Orro and Kisrer A. J. Ind. Eng. Chem. 1948 40 345. 

[5] Hupesranp J. H. Solubility of Non- Electrolytes 1936 Reinhold. 

E. G. and D. Ind. Eng. Chem. 1947 39 1329. 


(7] Barker J. A. Austral. J. Chem. 1953 6 207. 


Horsiey L. H. Ind. Eng. Chem. Analyt. 1947 19 508. 
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iso-octane- furfural benzene | measured | calculated 
4 zB rc mm Hg | mm Hg az = mole fraction in a liquid 

y = activity coefficient 

0-0798 0-7972 0-1230 69-2 69-5 Subscripts 

0-1285 0-6450 0-2265 17-1 79-5 .  B-— furfural 

0-8931 0-0545 0-0525 60-4 60-3 C= ‘“ C — benzene 

0-8077 0-0720 0-1203 72-0 7146 D-«= ‘. D — carbon tetrachloride 

0-6042 0-1439 0-2519 81:3 81-0 E= - E — acetic acid 

(1) 
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The equations of change and the macroscopic mass, momentum, 
and energy balances 


R. Byron Birp 
Department of Chemical Engineering, University of Wisconsin, Madison, Wisconsin 
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Abstract—It is well-known how the Bernoulli equation of mathematical hydrodynamics is 
obtained by forming the scalar product of the velocity vector with the Euler equation of motion 
for an inviscid fluid and integrating the resulting equation along a streamline (StREETER [8]). It 
is here shown how the “ engineering Bernoulli equation " — containing the mechanical work 
term and the friction loss term — may be obtained by a similar integration of the Navier-Stokes 
equation of motion for a viscous fluid over the volume of a flow system. For the sake of completeness 
it is first shown how the macroscopic mass, momentum, and energy balances used in engineering 
are obtained from the “ equations of change ™ of theoretical physics. In these derivations the 
VOL role played by the moving surfaces is taken into account explicitly, and the fluid flowing in the 
6 ' system may consist of several chemical species with chemical reaction occurring. The development 
956/57 leads to the unsteady state expressions for the macroscopic balances and the Bernoulli equation. 
The correct types of average velocities to be used in the macroscopic balances arise naturally 

in the course of the derivations. 


Résumé—L’obtention de l'équation de BeRNouLLI en mécanique des fluides est bien connue & 
partir du produit scalaire du vecteur vitesse par l’équation d’EULER du mouvement d'un fluide non 
visqueux que intégre le long d'une ligne de courant (Streeter 8). L’ auteur indique ici comment 
on peut obtenir par une intégration analogue de l'équation de Navirer-Sroxes, équation de 
mouvement d'un fluide visqueux, l'équation de utilisable pratiquement, contenant 
les termes de travail mécanique et de perte de charge par frottement. Dans un souci de précision 
l'auteur montre comment il obtient a partir des équations d’échange de la physique théorique, 
les bilans massiques, les bilans des moments et les bilans énergétiques. I| tient compte aussi 
dans ses calculs ; de fagon explicite, du réle joué par les surfaces en mouvement, et du fait que 
le fluide qui s‘écoule peut contenir plusieurs espéces chimiques en réaction. Le développement 
conduit & des expressions d’états de transition pour les bilans macroscopiques et l’équation 
de BrERNOUILLI. 

Les expressions correctes des vitesses moyennes dont on se sert dans les bilans macroscopiques 
apparaissent naturellement en cours de calculs. 


THe Equations or CHANGE FOR A which moves with the local fluid velocity (Birp 
MULTICOMPONENT SINGLE-Puase [1] and Hirscurevper et al. [8].4) 


Tue differential equations describing the motion 
of a fluid containing v chemical species are the 
“ equations of change which may be derived + 1% 

from non-equilibrium statistical mechanics 

(Hirscure.per et al. [3]). These equations in the i= 1, 2,3,...¥ (1)* 
forms given here are applicable to both laminar — ixplictt expressions for the various vector and tensor 
and turbulent flow, and apply to non-Newtonian operations in cartesian coordinates are given at the end 
as well as Newtonian fluids. They may be written of the artiele. 

in two equivalent forms: In the first form they * Addition of all » equations of the form of eq. (1) gives : 
indicate how the mass density of the ith chemical 


(la) 
species p;, the velocity v, and the total energy per Dt ° 


unit mass E change in a small element of fluid which is the overall equation of continuity. 
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In the second form, derivable from the first, the 
equations are written so as to show them as laws 
of conservation of mass, momentum, and energy 
applied to a small volume element (fixed in space) 
through which the fluid is flowing : 


Eq. of Continuity : 


=—(V {piv + t+ Tig 


Eq. of Motion : 


pv = —(V+{ovv + + pk, (5) 


Eq. of Energy : 


a A 
{pEv +(m+v) + q}) + pQ (6) 
+ Several alternate and completely equivalent formula- 
tions of the equation of energy balance are available ; one 
is in terms of internal energy and the other in terms of 
temperature : 


= ~(¥ + q) —(w: + py 


op 
+( vy) 


U; op 
+ ‘jp + val [m+ (3b) 
i=l 


in which the derivative )p/)T is taken at constant com- 
position and volume, and U; and V; are partial molal 
quantities. For single component, constant pressure 
systems eq. (3b) may be simplified to : 
a DT A 

which is the basis for most engineering calculations. The 
term (rt: Vv) is the heat produced by viscous dissipation 
(Lamp [4]) which may generally be neglected in flow 
problems. 


The physical significance of the “ microscopic 
balances” in eqs. (4), (5), and (6) is, term by 
term: The equation of continuity states that the 
local mass density of the ith species within a 
small region changes because of (i) the net influx 
of mass of i across the surface of the region by 
mass flow, (ii) the mass flux by diffusional pro- 
cesses (ordinary, thermal, pressure, and forced 
diffusion), and (iii) the production of the ith 
species by chemical reaction. The equation of 
motion states that the momentum within a small 
region changes because of (i) the flux of momentum 
across the surfaces by mass flow, (ii) the flux of 
momentum across the surfaces by pressure and 
viscous forces, and (iii) the external body forces 
(which are assumed to be the same for all species 
and derivable from a potential : F, =— Ve). 
The equation of energy states that the total energy 
within a small region changes because of (i) the 
flux of energy across the surface by mass flow, 
(ii) the work done by the fluid in the region on its 
surroundings (both pressure-volume work and 
work asscociated with viscous forces), (iii) the flux 
of energy entering the system by heat conduction, 
radiation, and mass transfer, and (iv) the heat 
produced within the fluid by nuclear radiations.* 
Equations (4), (5), and (6) are all of the form : 


(Accumulation) = (Input — Output) 
+ (Production) 


applied to a fixed differential volume element 
through which the fluid is flowing. 

The equations of change together with boundary 
and initial conditions, the thermal equation of 
state (p = p(p, T, p;)), and the caloric equation 
of state (U = U(p, T,p,)) give the complete 
description of the flow system. In order to caleu- 
late concentration, velocity, and temperature 
profiles, one has to substitute expressions for 
j;, t, and q in terms of the transport coefficients 
(coefficients of diffusion, viscosity, and thermal 
conductivity) and the gradients of the macroscopic 


* The energy release accompanying chemical reaction 
is not included as a separate contribution in eqs. (3) and 
(6), inasmuch as it is inherently included in the internal 
energy contribution to BE. The heat of reaction does 
appear explicitly when the energy equation is rewritten 
in terms of temperature (see r;, terms in eq. (3b)). 
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Eq. of Motion : 

Eq. of Energy : 

p _ — (wv) +0 (ot 
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The equations of change and the macroscopic mass, momentum, and energy balances 


properties (concentration, velocity, and tempera- 
ture). The equations thereby obtained are useful 
for laminar flow calculations ; their time-averaged 
counterparts are used as the starting point for 
turbulent flow calculations. 

The equations of change written in eqs. (4), (5), 
and (6) may be integrated over the volume of any 
macroscopic engineering flow system, such as the 
one shown in Fig. 1. Thereby one obtains the 
macroscopic (or “ overall ’’) mass, momentum, 
and energy balances used in the analysis of 
engineering systems. One may also take the 
scalar product of v with the equation of motion 
in the form of eq. (2) and integrate the resulting 
equation over the volume of a macroscopic system 
to obtain the Bernoulli equation in the form 
usually used in engineering. It is well known 
that these overall balances give no information 
concerning the concentration, velocity, and tem- 
perature gradients in a system but simply con- 
nections between properties of the inlet and outlet 
streams. 


Q, 


(a) (b) (c) 


Fic. 1. Macroscopic flow system with mechanisms for (a) 
the addition of heat Qy through fixed surface Sy and volume 
heat source Q¢ resulting from subatomic radiations, (b) 
transfer of mechanical work W to surroundings by means 
of moving surface S,,, and (c) both homogeneous chemical 
reactions Rj, and heterogeneous (catalytic) reactions Ry. 


In the integration of the equations of change 
use is made of the Gauss divergence theorem to 
change the integration over the volume (V) of 
the divergence of a vector to an integral of the 
normal component of the vector over the system 


surface (S). In doing this three types of surfaces 
have to be distinguished : the fixed wetted solid 
surface in the system (S,), the movable wetted 
surface which allows for the transfer of energy 
between the macroscopic system and its sur- 
roundings in the form of mechanical work (S,,), 
and the open end surfaces at the control surfaces 
(S,, S,). These different surfaces have to be 
treated individually. Throughout it is assumed 
that the fluid properties (e.g. density) do not vary 
appreciably over the control surfaces S, and S,. 
The control surfaces are chosen in such a way 
that the flow velocity at S, and S, is, aside from 
turbulent fluctuations, parallel to the pipe axis 
and normal to the control surfaces. It is further 
assumed that transport of mass, momentum, and 
energy across S, and S, by molecular processes 
may be neglected in comparison with the transport 
by mass flow. 


2. Macroscopic Mass BALANceE 


Integration of eq. (4) over the volume of the 
macropscoic flow system V gives :* 


ay — 
dt 


[ov frgav (7) 
Vv 


The first two integrals on the right hand side 
may be transformed into surface integrals by 
means of the divergence theorem. The integral 
on the left hand side may be rewritten by taking 
the time derivative outside the integral sign. In 
so doing, however, one must remember that 
because of the moving parts within the flow 
system, the shape of the volume V is changing 
with time, and hence a term must be subtracted 
which accounts for the local gain or loss of mass 
by virtue of the motion of the surface. One there- 
by obtains : 


* The mass flux j; is the local flux of component i in 
gms cm sec~! referred to a co-ordinate system which is 
moving with the local mass average velocity v. The 
relation of this flux to other fluxes (including the flux N; 
more usual in chemical engineering) is given in ref. [1), 
Table III. 
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Here v,, is the velocity of the fluid normal to the 
surface of the system, and »,,,, is the component 
of the velocity of the moving surface which is 
normal to the surface (see Fig. 2). Hence eq. (8) 
becomes : 


Fic. 2. A differential element of surface dS of the moving 

surface of the macroscopic system moves with an instant- 

aneous velocity v,,; the component of v,, normal to dS 
is 


M, Jo — Umn) aS 


8, +8, Sy Vv 
= 5, — p25.) 
| Jin dS + Ry + Ri 
Si 
In the first term on the right hand side use has 
been made of the assumption that p; is constant 


over the cross-section S, and S,. The integral of 
Jin, over S, and S, may be neglected in comparison 


> (9)t 


+ In what follows. we use 0, 0%, for averages over a 
cross-sectional area ; specifically for a circular tube these 
are defined by 

2a 2e R 


where n is 1, 2, or 3. 
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with the first term which represents mass trans- 
port of the ith species by the overall flow motion. 


The integrals Ry and Ry are the net rates of 
production of the ith species by surface reactions 
on S, and by homogeneous reactions in the fluid 
phase respectively. If now R; = Ry + Ry and 
J a = PadyS), 


M,=-—AJ,+R, i=1,2,...» (10) 


where AJ; = Jig — Jy. Addition of all the v 
equations represented by eq. (10) gives : 
M=—AJ (11) 


Equations (10) and (11) are the familiar macro- 
scopic mass balances for individual chemical 
species and for the fluid as a whole. For steady 


state operations M,; and M are zero. 


8. Macroscopic MomENTUM BALANCE 
Integration of eq. (5) over the volume gives :* 


| 


> (12) 
(V- nav 


The left hand side may be transformed as des- 
cribed in Section 2. And the first three volume 
integrals on the right hand side may be trans- 
formed into surface integrals by means of modi- 
fications of the divergence theorem (MILNE- 
Turomson [5], Morse and Fesupacn [6)). 


* The element 7,, of the second-order pressure tensor 
w is the force per unit area acting in the 2-direction on an 
elemental fluid surface perpendicular to the y-direction. 
The diagonal elements 7,,, tyy, correspond to normal 
forces, and the non-diagonal elements "ay = "ye 
= and = correspond to shearing forces. 
The diagonal elements include both hydrostatic pressure 
forces and viscous forces: = + Tg,; the non- 
diagonal elements include only viscous forces : 7, = Tay: 
This is summarized by the tensor relation w = p& + 1, 
in which 8 is the unit tensor (5,, = 5,, = 3,, = 1, all 
other components are zero). In eq. (12) w is split up into 
p and ¢ and use is made of the fact that (Vi: w) = Vp + 
7). 
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Here n is the normal unit vector directed out- 
wardly from the system surface. The left hand 
side is the time rate of change of the total 
momentum P of the fluid in the system. The 
first term on the right hand side vanishes on all 
solid surfaces; the contribution at S, and S, 
represents the rate of flow of momentum into 
and out of the system by virtue of fluid flow 
alone. The second term on the right gives a 
contribution to the influx and efflux of momentum 
due to pressure forces at S, and S,, and also a 
contribution to the drag force on the equipment 
on S, and S,,. The third term also contributes 
to the drag force on the equipment due to viscous 
forces on the solid surfaces S,and S,,; in addition 
this term includes momentum transport into and 
out of the system at S, and S, by viscous action — 
a contribution almost always safely neglected. The 
fourth term gives the total force exerted on the 
fluid in the system by external body forces such 
as gravity. Hence one obtains finally the macro- 
scopic momentum balance : 


P = — A(pv® S + pS) — Pag + Fo (14) 


or in terms of the mass rate of flow : 


> (13) 


P = -a(5 J + pS) — Fi + Fin (15) 
Here S indicates a vector whose magnitude is the 
surface area of the cross section and whose direc- 


tion is normal to the cross-section (chosen to point 
in the direction of flow!). The vector J is the 
mass rate of flow vector similarly defined. The 
negative of F,,,, is just the “reaction force ” 
offered by the equipment. Note that the average 
value of the square of the velocity appears 
properly in the momentum balance (Rouse and 
Howe [7)). 


4. Macroscopic ENerGy BALANcE 
Integration of eq. (6) over V gives : 


Jov-wav + 

The first three volume integrals on the right hand 
side may be transformed into surface integrals as 
before. 


The left hand side is the rate of change of the 
total energy in the system. On the right hand 
side the first integral is zero except on S, and S, 
and represents the energy flowing in and out with 
the fluid. The second integral vanishes on S,; its 
value on S,, gives the rate at which work is done 


on the external surroundings, W, resulting from 
both pressure and viscous forces on the moving 
parts: on S, and S, it contains a term pv, (the 
pressure-volume work) and a term (r+ v), (the 
shear work), the latter being generally neglected. 
The third integral gives the heat added to the 


system over the solid surface Qi there is also a 
contribution at S, and S, which represents the 
heat input and outflow by conduction in the fluid 
stream, this contribution generally being quite 
small in comparison with the heat transport by 
bulk flow. Finally the fourth integral represents 
the volume heat production Q,. The surface 
integral over E can be shown as being made up 
of internal energy U, kinetic energy }v%, and 
potential energy @ (all on a per unit mass basis). 
Then eq. (17) finally becomes : 


E = — A + + po OS + pos) 
+Q-W 


«as) 


+ 
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in which E = U + K + }and Q = Q,+ Q,. Or, 
in terms of the mass rate of flow, the macroscopic 
energy balance is : 


E=—A(U + pP + + 6) J + Q—W 


Here it is seen that the average value of the cube 
of the velocity enters quite naturally (Rouse and 
Howe [7)). 


5. Tue “ENGINEERING BERNOULLI 
EQuATION” 
The formation of the scalar product of v with 
eq. (2) gives (StreeTeR [8]) after some manipu- 
lation : 


= —(v+ Vp) 
—(v+ (V+) 


} (20) 
This is an energy equation involving only mech- 
anical terms —that is, heat and internal energy 
do not appear. The left hand side may be re- 
written by using the equation of continuity to 
give terms (a) and (b) below; and the second 
term on the right hand side may be split up into 
two terms (d) and (e), provided that use is made 
of the fact that + is a symmetric tensor. Thereby 
one obtains : 


(dpc?) = — v) — (v + Vp) 
(a) (b) {e) 
—(V+(r+v)) + Vv) + p(v-F,) 
(@) (e) 


* For steady-state J, = J, and eq. (19) can be written 
on a per unit mass basis by dividing through by J. 


A(U + ph + + 6) 


It should be noted that Q does not include the heat pro- 
duced by chemical reaction inasmuch as that effect is 
automatically included in the calculation of the internal 
energies U, and Us of the inlet and outlet streams. 


(19a) 


This equation is next integrated over V to give : 


Term (a): 
d 
J 


Here K is the rate of change of the total 
kinetic energy within the system 


Term (b) : 


J 


— 
Sm 


Sm d 


The integral over S,, just cancels the same 
integral in eq. (22); the other term gives the 
input of kinetic energy at surface 1 minus the 
output at surface 2. 


Term (d) : 


— [(v-(r-w) av= -| (r+Vv),dS = — W, (24) 
Vv Sm 


The contributions of the surface integral over S, 
anc S, are neglected as was done in going from 
eq. (17) to eq. (18); W, is the rate of transferring 
mechanical energy to the surroundings because 
of the viscous forces. 
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Term (e) : 


+ &, (25) 


This is the total rate of energy loss by viscous 
dissipation within the system (Lamp [4]) — that 
is, the friction losses resulting from the degrada- 
tion of mechanical energy into heat. 


Term (f) may be evaluated by introducing the 
potential energy, which is related to 


[ — | ove vdav 
Vv Vv 


The first of these integrals may be transformed 
by using the divergence theorem and the second 
may be rewritten by using the equation of 
continuity (and the fact that @ is time- 
independent) : 


The first term is just the influx and efflux of 


potential energy, and @ is the rate of change of 
the potential energy of the entire system. 


Term (c) requires special consideration. If the 
fluid is of constant temperature and composition, 
then use can be made of the fact that (1/p)dp=dG 
and the development parallels the treatment of 
term (f): 


Pp 


d a 
Ss 


A(pGas) — A4-W, 


The first term gives the rate of influx and efflux 
of G; the second term is the time rate of change 
of the Helmholtz free energy A = G — pV for 
the whole system ; and the third term is the work 
done by pressure forces on the moving parts. 

Combination of the results of eqs. (21) through 
(27) gives then finally : 


(28) 


where W = W, + W,. This is the unsteady state 
engineering Bernoulli equation for an isothermal 
system of constant chemical composition. 

For an isentropic system of constant composi- 
tion, term (c) is evaluated by noting that 
(1/p)dp = dH. A similar procedure leads to the 
result 


K+640 = — (29) 
For a system which is neither isothermal nor 


isentropic or in which chemical reactions are 
taking place no such simple results may be 
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| ovdas — [ av 
dt 
d A 
A GdS — — V 
F, by F,= — VO: a | (27) 
+ a8 + Pav 
s v 
VOL. = | (0, — 
6 s 
— 
s 
— [w,das 
§ 
+ 
(26) 
s 
= — — 
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written down. The two cases described by eqs. 
(28) and (29) are those usually considered in 
engineering text (Coutson and Ricuarpson [2)). 

For steady state operation J, = J, = J, and 
the engineering Bernoulli equation may be 
written on a per unit mass basis by dividing 
through by J to obtain : 


a 


a Po 
aye +ad+ + W+BE,=0 (80) 
Pi 


For the isothermal, eonstant-composition process 
the line integral reduces to AG ; for the isentropic, 
constant-composition process the line integral is 
AH. For any other process the best one can do 
is to evaluate the line integral over the succession 
of thermodynamic states encountered in the flow 


A = Helmholtz free energy 
C, = heat capacity at constant pressure 
C, = heat capacity at constant volume 
E = total energy 
E,, = energy loss due to friction 
F, = volume force = — V® 
F.., = total external force 
F = total drag force 
=A+pV 
H = enthalpy 
j; = mass flux of ith component with respect to the 
mass average flow velocity (Zj; = 0) 
J = mass rate of flow 
K = total kinetic energy of macroscopic system 
M = total mass of fluid in macroscopic system 
M, = mass of the ith species in the macroscopic flow 


p = static pressure 
P = total momentum of system 
q = heat flux vector 
q, = normal component of q 
Q = heat transferred across surface of macroscopic 
system 
Q, = volume heat source in macroscopic system 
Tig = mass of ith chemical species produced by homo- 
fe eneous chemical reaction (Zr; = 0) 
R, = total rate of production of mass of the ith species 
in the macroscopic system 
S = surface of flow system 
Sy. Sm = fixed and moving surfaces of the system 
S,, Sg = cross-sectional area of flow system at 1, 2 
t = time 


R. Byron Birp 


NOTATION 


from S, to S,, the integration being performed 
over some representative streamline in the 
system. This implies that one know, for example, 
how p varies with p, 7, and composition (specifi- 
cation of the equation of state) and how T and 
the composition vary with p (specification of the 
path of the line integral). 
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U = internal energy 
Vv = velocity 
v, = normal component of v 
V = volume of flow system 
o = velocity averaged over a cross-section 
Vv, = Velocity of surface of a moving part 
Umn = normal component of v,, 
W = work done by system on its surroundings 
W, = work resulting from pressure forces, viscous 
forces on moving parts of the equipment 
4 = quantity evaluated at S, minus quantity evalu- 
ated at S, 
8 = unit tensor 
7 = pS + 7 = pressure tensor 
p = mass density = (1/7) = Zp, 
p; = mass density of ith chemical species 
+ = shear stress tensor 
@® = potential energy 
Marks above symbols : 
«a = quantity per unit mass of fluid 
— = quantity averaged over cross-section of flow 
system 
+ = time rate of change of quantity, pertaining to 
the macroscopic system 
Subscripts : 
f = quantity associated with fixed surface 
m = quantity associated with moving surface 
¢ = quantity associated with the fluid 
n = normal component of a vector 
v = quantity associated with viscous effects 
Pp = quantity associated with pressure effects 
i = ith chemical species 


vol 
6 

system 
n = outwardly directed normal unit vector 
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Vector and Tensor Operations in Cartesian Co-ordinates* & = second order unit tensor with three compo- 
Scalars are indicated by light face symbols nents 3;; = 1 and six components 3, = 0 
Vectors are indicated by bold face Roman symbols 

Tensors are indicated by bold face Greek symbols + the gradient of 
i,j,k = unit vectors in z, y, z directions de aw _ 
v = iv, + ju, + kv, = vector with 3 components 
Up 
7 = second order tensor with 9 components +,,, es ee 
Tay Tax ete.t dz oy dz 
vv = dyadic product with 9 components v,0,, 0,0,, > > > 
ete. (V+ pvV), = — + — + — pu, 
az dy dz 
v? = (v-v) = 0,3 + + 


(v+ Ww) = + ty, + = scalar product of two (v + =(v: Up) = 0, +0, +0, 
vectors dz dy dz 


(T+ = + TayYy + Vv) =z rap dp) (a, B = 2, y, 2) 
a 
* See ref. [1] for a tabulation of the more difficult Vv = dyadic product with nine components dv,/dz, 
operations in cyclindrical and spherical co-ordinates. duy/dz, dv,/ dx, dv,/dy, ete. 
VOL. + All tensors used in this article are “symmetric” — that + V)p = the “ sul tial derivative 
Ti = 


of p” 
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Abstract—An improved form of the Bearrie volume-explicit equation of state for carbon 


dioxide is described, containing two additional terms. Tables of compressibility, fugacity, and 
water-solubility of carbon dioxide are presented for the region 0-36 atm. and 0-100°C. 


Résumé—L auteur présente une forme améliorée de l'équation explicite d'état de Bearrie, 
équation d'étar du gaz carbonique avec deux termes additionnels. Il donne des tables de com- 
pressibilité, fugacité et solubilité de l'eau de anhydride carbonique pour les régions de 0-36 atm. 


et 0-100°C, 


INTRODUCTION 


For the design of a carbon dioxide absorber 
operating at pressures up to 36 atm. [1] it 
has been found necessary to revise existing 
data on the compressibility, fugacity, and water- 
solubility of carbon dioxide for pressures in the 
range 0-36 atm. and temperatures in the range 
0-100°C. It is hoped that the revised data in this 
paper will prove useful to other workers. 

To calculate the fugacity it was necessary to 
obtain an equation of state, explicit in volume. 
The only volume-explicit equation available in 
the literature for the required range of tempera- 
ture and pressure was the Beattie equation [2], 
which had been obtained by inversion of the 
Beattie-BripGeMAN pressure-explicit equation 
[3]. Both equations were based upon the data of 
Anprews [4] for the compressibility of carbon 
dioxide. When the Beartie equation was tested 
by using the more recent compressibility data of 
MicHeELs and Micue [5], the average deviation 
of the calculated from the experimental volumes 
was 2.4%, much too great an error for the 
accurate computations now required. 

ZeELVINSKI [6] determined the water-solubility 
of carbon dioxide at certain temperatures in the 
range 0-100°C and pressures up to 90 atm., while 
Wiese and Gappy [7] gave the solubilities of 
carbon dioxide for a single pressure (25 atm.) 
within the range now required, and at seven 
temperatures between 18° and 100°C, There are, 
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in addition, the solubility data (1 atm.) of Bour 
and Bock [8], quoted in most handbooks, Since 
all these data are expressed in different units, it 
was decided (i) to evaluate their accuracy and 
compile a table of solubilities suitable for design 
purposes, and (ii) to test the validity of Henry’s 
law (and its fugacity correction) so that the 
chemical engineer might have a guide to the use 
of this law for the system carbon dioxide-water. 


EQUATION OF STATE FOR 
CarBON DIOXIDE 

Anprews [4] obtained his data by comparing 
the volumes of carbon dioxide and air at the 
same temperature and pressure. The com- 
pressibility of air was unknown in 1876, and 
ANDREWS assumed it to be unity; however, it 
was possible to refine his data by applying a 
correction for the compressibility of air, using the 
data of HoLporn and Orro [9]. Based upon the 
corrected ANDREws’ data, the volumes calculated 
from the Beattie equation had an average 
deviation of 1-6%, with a maximum of 6.6%. 

A review of equations of state shows that 
volume-explicit equations require more empirical 
constants than pressure-explicit equations to 
achieve the same accuracy. The Beattie equation 
has therefore been modified, by the addition of 
two further terms, to fit the more recent experi- 
mental data [5]. It can be expanded into the 
following volume-explicit form : 
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RT 61-0102 660,000 ) 
P 4-35126  69141-6 
1 
aT 0.007579 + (1) 
RT T® 


When the most recent compressibility data for 
carbon dioxide [5] are used to test equation (1) 
the calculated molal volumes show an average 
deviation of 2-4°%, and a maximum deviation of 
11-7%, from the experimental values. Further, 
the deviations are such that the calculated 
volumes are always larger than the experimental. 
Graphical and analytical examination of these 
deviations has now shown that the Beatriz 
equation can be adjusted to fit the Micuets and 
Micue.s data by the addition of two further 
empirical terms. The deviations are found to be 
a function of both temperature and pressure ; 
hence, one term, — 2-47 x 1077/7", is now 
added to allow for the temperature effect, and a 
second term — (P/ RT’)? (2-69354 10" to 
allow for the combined effects of temperature and 
pressure. This leads to the following volume- 
explicit equation : 


RT 61-0102 
V = 0.10476 + 
660,000 2-47 x 
(2) 
P 435126 69141-6 
P 5Q02-39 2.69854 x 10'° 
RT Ts J 


Equation (2) can be converted into the com- 
pressibility form by multiplying through by P/ RT 
to obtain : 


PV 
RT 
61-0102 660,000 2-47 10?" 

[o-r0476.- oe 

P \? 4-35126 69141-6 
[0.007579 

P 5002-39 2.69354 x 10° 
TS 
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These equations of state can be applied to 
carbon dioxide, over the pressure range 0-36 atm. 
and the temperature range 3-100°C, with an 
average deviation of 0.17% and a maximum 
deviation of 0.35%. They can also be applied, 
with the same deviations, over the temperature 
range 0-8°C and the pressure range 0—33-4 atm. 
The vapour pressure of liquid carbon dioxide at 
0°C is 33-4 atm., so that within the narrow range 
0-3°C the substance will be liquid at pressures 
above 33-4 atm. Equations (2) and (3) have been 
derived solely for the purpose of interpolation so 
that accurate compressibility coefficients could 
be tabulated. These coefficients are given in the 
Table at equal temperature and pressure incre- 
ments. 

It should be noted that equations (2) and (3) 
also agree well with the older ANDREws data [4], 
with an average deviation of 0.38% and a 
maximum deviation of 0-53%. However, the 
deviations are always positive, indicating that 
there is some constant error in the ANDREWs 
data which makes the volumes consistently 
smaller than those reported by Micnets and 
MICHELS. 


COMPRESSIBILITY 


The compressibility of carbon dioxide, as calcu- 
lated from equation (3), is given in the Table for 
pressure increments of 2 atm. up to 36 atm., and 
for temperatures of 0°, 10°, 15°, 20°, 25°, 85°, 50°, 
75° and 100°C. For design purposes, intermediate 
values can be interpolated; alternatively, a 
compressibility chart can be constructed. 


FuGAcItTY 


The fugacity coefficient of carbon dioxide is 
obtained by substituting equation (2) into the 
following integral : 


V 
laf = | —— 
nf RT dP (4) 


The resulting definite integral is integrated term 
by term to obtain the following expression for 
the fugacity coefficient : 


G. Hovcuron, A. M. McLean, and P. D. Rricuie 


0.10476 — °1:0102 The fugacity coefficients (f{/P) calculated from 
P RT| a. equation (5) are summarised in the Table. 


660,000 2-47 x 1 
TR oF Carson Dioxipe 


P\? Before the results of Ze.vinsku [6], Wiese and 
(er) Gappy [7], and Bour and Bock [8] can be com- 
0.0087895 — 2.17563 4 34570-8 pared, all the data must be converted into the 

T T? same units of solubility. The units used here are 


+ ( P x the Bunsen absorption coefficient («), and the 


RT T? T mole fraction of dissolved carbon dioxide (7). 


Table 1. Compressibility, Fugacity, and Water-Solubility of Carbon Diovide 
Z = PV/RT = compressibility coefficient ; {/P = fugacity coefficient ; 


a = Bunsen absorption coefficient, volume of gas calculated at N.T.P. dissolved by 1 volume of water at 
a certain temperature and partial pressure P of CO, ; 


2 = mol. fraction of dissolved CO, in water at a certain temperature and partial pressure P of COg,. 


oc 10°C 


j/P | 
0-99263 0-99266 0-99362 
0-98527 098536 P O-98724 | 
0-97055 0-97098 097470 
0-95576 0-95676 . | 096230 
0-94082 094275 0-95004 
0-92567 0-92886 0-93788 
0-91022 O-91512 0-92584 
0-89440 090146 0-91388 
O-87814 0 88792 0-90198 
0-86136 0-87442 0-89016 
0-8-4400 086097 0-87836 
0-82596 084753 086657 
080719 0-83413 085483 
O-78760 082073 084305 
76712 080731 083130 
074567 0-79387 0-81950 
72319 0- 78039 080770 
0-79583 
_ 0-78393 


) 


P 
atm. 

1 

2 

4 

6 

8 
10 
12 
14 
16 
18 
20 
22 
24 
26 
28 
30 
32 
34 
36 


20°C 


0-99400 0-99400 099437 099436 
0-98800 0-98806 0-98874 0-98880 
0-97602 0-97630 0-97750 O-97772 
0-96401 096464 0-96624 0-96678 
0-95193 0-95314 0-95491 0-95596 
0-93971 O-04174 094548 094524 
0-92732 0-93042 0-93190 0-93458 
0-91470 0-91918 . 0-92013 0-92402 
0-90181 0-90796 . 0-90812 0-91350 
0-88860 0-89690 O-89584 0-90302 
O-87501 0-88580 O-88324 0-89258 
0-86101 O-87474 P O-87029 0-88218 
0-84654 0-86370 0-85693 0-87178 
083156 0-85263 0-84312 0-86140 
O-81601 0-84160 082883 0-85100 
0-79985 0-83053 081401 0-84058 
0-78303 0-81941 0-79862 0-83013 
0-76549 0-80830 0-78262 0-81967 
0-74720 O-79713 0-76596 0-80917 
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1218 (0-085 
2-41 | 1-046 

471 3-80 

691 5:56 

9-00 7-23 

10-98 8-81 

12-86 10-30 

14-63 11-70 

16-29 13-02 

17-85 14-24 

19-30 15-38 

20-64 16-44 

21-88 17-41 

23-01 18-29 

24-04 19-09 

24-96 19-80 

25-77 20-44 

26-48 20-98 

| 27-08 21-45 

15°C 

0-855 0-692 
1-697 1-374 

3-35 2-70 

494 | 3-99 

6-50 5-24 

8-00 | 6-44 

9-45 7-60 

10-85 | 8-72 

1221 |) OBI 

13-52 10-84 

14-78 11-84 

15-99 12-80 

17-138 

18-26 | 14-59 

1932 «15-43 

20-34 16-22 

21-30 16-98 

22-22 17-70 

22-99 18-30 

| | 
184 
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25°C 


Fr 


0-81259 
0-79793 
0-78270 


0-99472 0-99532 
0-98946 “ 0-99062 
0-97904 098136 
0-96874 0-97218 
0-95854 096310 
0-94844 0-95410 
0-92842 093624 
0-91852 092738 
0-90864 0-91856 
0-89882 0-90980 
088900 0-90106 
0-87918 0-89230 
0-86940 0-88356 
0-85960 O-87484 
084976 0-86610 
083995 0°85735 
083010 0-84857 
082021 0-83979 


50° 75°C 


1 

2 

4 

6 

8 
10 
12 
14 
16 
18 
20 
22 
24 
26 
28 
30 
82 
34 
36 


0-99603 
0-99206 
0-98414 
0-97623 
0-96831 
0-96035 
095234 
0-94426 
0-93608 
0-92779 
0-91938 
0-91081 
0-90208 
0-89316 
0°88403 
0-87467 
086507 
085520 
0°84505 


0-99602 099696 
0-99208 . 0-99393 99394 
0-98424 . 0-98788  0-98794 
0-97646 . 0-98184 0-98196 
096876 . 0-97581 0-97604 
0-96112 0-96976 0-97014 
0-95352 0-96369 0-96432 
0-94596 . . 0-95759 0-95850 
0-93844 . 0-95145 0-95268 
0-93098 0-94525 0-94692 
0-92350 . 0-93899 0-94116 
0-91606 . 0-93266 0-93542 
090866 . 0-92625 0-92972 
0-90124 . 0-91973 0-92402 
0-89384 . 0-91312 0-91832 
0-88644 . 0-90638 0-91260 
0-87902 . 0-89952 0-90688 
0-87160 . 0-89253 0-90122 
0-86417 0-88538 0-89550 


100°C 


0-99766 
0-99532 
0-99066 
0-98601 
0-98137 
0-97678 
0-97207 
0-96740 
0-96272 
0-95801 
095327 
0-94849 
0-94367 
0-93880 
0-93388 
0-92889 
0-92383 
0-91871 
0-91350 


| 35°C 
Z | | 2 x 108 Z | f/P | x 108 
0-99470 0-581 0-473 
0-98941 1-159 0-943 
0-97885 2-30 1-868 
096827 3-42 2-78 
0-95764 452 3-67 
094691 5-60 4-54 
0-93607 6-66 5-39 
0-92506 7-70 6-23 
091385 8-72 7-05 
090240 9-72 7-85 
0-89069 10-70 8-64 
087866 11-66 9-40 
086629 12-60 10-15 
085353 13-52 10-89 
084035 14-42 11-60 
0-82672 15-30 12-30 
16-16 12-99 
| 17-00 13-65 
17-82 14-30 
956/57 
9-299 0-248 
| 0-597 0-495 
| 1-187 0-984 
1-768 1-465 
2-34 1-941 
2-91 2-41 
| 3-47 2-87 
4-02 3-33 
4:56 3-78 
5-10 4-22 
5-63 4-65 
6-15 5-08 
6-67 5-50 
717 5-92 
7-67 6-33 
8-16 6-73 
8-65 7-13 
9-12 7-52 
| 9-59 7-90 
099764 0-221 0-187 
099534 0-442 0-373 
0-99070 0-881 0-743 
098606 1-317 1-111 
098150 1-751 1-477 
097692 2-18 1-841 
0-97240 2-61 2-20 
0-96786 3-04 2-56 
096338 3-46 2-92 
0-95890 3-89 3-27 
0-95442 431 3-62 
094996 4-72 3-97 
094554 5-13 4:32 
0-94108 5°55 4-66 
0-93668 5-96 5-01 
0-93224 6-36 5-35 
0-92784 6-77 5-68 
0-92342 7:17 6-02 
091898 7-57 6-35 
135 


G. Hovucuron, A. M. McLean, and P. D. Rircnie 


ZELVINSKII expressed his data at each tempera- 
ture by the equation : 


S=aP — bP* (6) 


(where S is the solubility in ml. of CO, calculated 
at N.T.P. dissolved in a quantity of water whose 
volume is 1 ml. at 0°C, and a and b are constants). 
Since ZELVINsk1! gives only a limited number of 
experimental solubilities at the temperatures 0°, 
12.48°, 25°, 50°, 75° and 100°C, solubilities at 
intermediate temperatures are obtained by 
graphical interpolation of the constants a and b. 
The accuracy of this procedure, checked against 
known experimental values, is found to be good 
within 2%. 

ZELVINSKI did not measure absorption co- 
efficients at 1 atm.; but these can be obtained by 
extrapolation of his data at higher pressures, 
using equation (6). The coefficients obtained by 
extrapolating to 1 atm. agree with those obtained 
experimentally by Bour and Bock, with an 
average deviation of 2-7% and a maximum of 
4-6%. 

The absorption coefficients at 25 atm. can be 
interpolated from ZeLviINsKil’s equations, and 
compared with the experimental values of WreBE 
and Gappy after correcting their data for the 
vapour pressure of water. The agreement is 
within an average deviation of 1-9% and a 
maximum of 5-7%. 

The accuracy of interpolation from ZeLvin- 
sku!’s equations has been checked by comparing 
interpolated values from equation (6) with actual 
experimental solubilities determined by him for 
various pressures [6]. The average deviation was 
1%, and the maximum 2.5%. 

The data given in the Table for the water- 
solubility of carbon dioxide have been obtained 
by interpolation of Ze.vinsku’s data and a 
detailed comparison of the results with the data 
of Wiese and Gappy and of Bour and Bock. It 
is concluded that the accuracy of these data is 
better than + 3%, and in agreement with other 
published data within the same limits. 


APPLICATION OF Henry’s Law 


Since the Bonr and Bock data for the absorption 
coefficients of carbon dioxide in water at 1 atm. 


are quoted in most textbooks, it was decided to 
apply Henry’s law to these figures and compare 
the results with values interpolated from the 
ZELVINSKII equations. In view of the maximum 
deviation deduced in the previous section (better 
than + 8%) it appears that deviations of more 
than 6% between calculated and experimental 
solubilities at the higher pressures can be con- 
sidered as deviations for Henry’s law. It is 
evident from the figure that calculations based 
on this law lead to solubilities consistently higher 
than the experimental values, the difference 
increasing as the pressure is raised and the 
temperature lowered. If 6%, is considered as the 
combined accuracy of the experimental data, 
Henry’s law can be applied only for pressures 
below ca. 5atm. within the temperature range 
0-100°C, 

To ascertain whether the non-ideal behaviour 
of the gas was responsible for the deviations from 
Henry's law, fugacity was used instead of 
pressure, again taking the Bour and Bock 
experimental data as a basis. The results of the 
fugacity correction are included in the figure. The 
effect of using fugacity instead of pressure is to 
correct Henry's law partially, and to produce 


Mole fraction (x) 


Fic. 1. Comparison of experimental solubilities of carbon 
dioxide in water with Henry's law and the fugacity cor- 
rection to Henry's law. 

———-——-_ Experimental solubility (Ze_vinsk11). 
Henry's law based on the absorption co- 
efficients of Bonur and Bock. 

—— — — Fugacity correction to Henry's law based 
on the absorption coefficients of Bour and 
Bock. 
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better agreement with the experimental data. At 
temperatures between 25° and 100°C, and at all 
pressures up to 36 atm., the fugacity correction 
can be applied to the Bonr and Bock data with 
an accuracy of better than 6%. At pressures 
below 10 atm., and at all temperatures between 
0° and 100°C, the fugacity correction gives an 
accuracy better than 6%. 

The conclusion reached is that below 25°C and 
above 10 atm. the system carbon dioxide—water 
is non-ideal. The deviations cannot be wholly 
attributed to departure from the ideal gas law, 
since the fugacity correction to HeNry’s law still 
shows considerable deviations, particularly at low 
temperatures and high pressures. It must therefore 
be assumed that the. activity coefficients of 
dissolved carbon dioxide differ considerably from 
unity. Evidence [10] from the rate of reaction of 
dissolved carbon dioxide with alkali indicates that 
less than 1%, of it is in the form of carbonic acid, 
which can therefore contribute little to the 
deviation of the solution from ideality. Further- 
more, temperature affects solubility much more 
strongly for carbon dioxide than for much less 
soluble gases such as hydrogen, oxygen, and 
nitrogen, which obey Henry's law over a wider 


range of pressure. This high temperature-coeffi- 
cient for the decrease in solubility of carbon 
dioxide shows that the heat of solution is large 
(relative to that of hydrogen, oxygen, or nitrogen) 
indicating compound formation between carbon 
dioxide and water. Carbonic acid being present 
only in small amount, the nature of such a com- 
pound is for the present a matter of conjecture ; 
it may be, for example, that hydrogen bonding 
between molecular carbon dioxide and water is 
responsible for the non-ideality of the solutions. 


Acknowledgment—The authors thank the 


Admiralty for permission to publish this paper. 


NOTATION 


V = volume of one g-mole of gas ; litres/g-mole. 

P = partial pressure of gas ; atm. 

T = absolute temperature ; °K. 

R = universal gas constant ; 0-08206 litre, atm./g-mole, 
°K. 

Z = compressibility of gas (PV /RT) ; dimensionless. 

f = fugacity ; atm. 

«x = Bunsen absorption coefficient ; volume of gas re- 
duced to N.T.P. dissolved by one volume of water. 

x = mol. fraction of dissolved carbon dioxide. 

S = solubility in ml of gas (reduced to N.T.P.) dissolved 
in a quantity of water whose volume is 1 ml at 0°C. 

a,b = empirical constants. 
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Flooding velocities in packed columns operating at reduced pressures 


H. SawistowskI 
Department of Chemical Engineering, Imperial College of Science and Technology, London, S.W.7 


( Received 30 October 1956) 


Abstract—The generalized graphical correlation for flooding rates in packed columns presented 
by Loso et al. [2] can be expressed in the form of an empirical equation as 


vere 


On the basis of this equation it has been shown that in a column operating at reduced pressures 


flooding will start at the top of the column. 


Résumé—L ‘auteur montre que le calcul des vitesses d’engorgement dans une tour 4 garnissage 
par la méthode graphique de Loso et al. peut s‘exprimer par nue équation empirique de la forme : 


Il a montré, & partir de cette équation, que dans une tour a pression réduite, engorgement 


INTRODUCTION 


Tue usual procedure for determining the flooding 
velocity in a packed column [1], [3] is to use the 
graphical correlation prepared by Lono, Frienp, 
HasuMa.t and Zenz [2], which is based on earlier 
work by Suerwoop, SuipLey and [4}. 
In this correlation a log-log scale is used and the 


w 


&* pt & pt pa 
plotted against (22) » where 
G 
Uy, = flooding velocity of the gas phase based 
on total column cross-section, 
G, = flooding rate expressed as mass flow of 
the gas phase per unit area of column, 
a = surface area of packing per unit volume 
« of column, 
g = acceleration due to gravity, 
« = void fraction of the packing, 
py = density of the liquid, 
Pg = density of the gas, 
zy, = Viscosity of the liquid, 
tw = Viscosity of water at 20°C (approx. 1 cen- 
tipoise), 
G = mass rate of flow of the gas phase, 
L = mass rate of flow of the liquid phase. 


group 


commencera a se produire au sommet de la colonne. 


In the case of a column operating at reduced 
pressure the problem is more complicated when 
the total pressure and the pressure drop are of a 
similar order of magnitude. The large variations 
of gas density along the column mean that it is 
not immediately obvious from the graphical 
correlation whereabouts in the column flooding 
will commence. The present paper is an attempt 
to solve this problem analytically. 


\ 


ht 


0-01 0 


+0 


Generalized correlation for flooding velocities in 
packed columns [2). 
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GENERAL EXPRESSION FOR FLOODING 
Rates iN Packep CoLumMNs 


The first step is to find an algebraic expression 
which fits Lopo’s graphical correlation. The shape 
of the curve (Fig. 1) indicates that substitution 
of a linear scale for logarithmic scale on the 
abscissa axis might transform the curve into a 
straight line whose equation would then be of 


& pt pc \uw 
For a given column, operating at a constant 
temperature, we can thus write 
=C 
| 


where A = —° (w=) = const. 
pr 


(2) 


If the experimental data are plotted as 


log (log versus le a straight 
line of slope n will be obtained if equation (1) 


is applicable, In Fig, 2, it is seen that points 
which are taken directly from Lopo’s curve do 
in fact lie on a straight line so that the values of 
C and n can be determined. The flooding rates 
are, therefore, correlated by the following 
equation : 


exe] 


Fic. 2. Evaluation of a general expression for flooding rates 
in packed columns. 


On substituting the value of A, the final expression 
Ga 


The accuracy of this equation is the same as that 
of the graphical correlation. 


In 


DETERMINATION OF THE POSITION OF 
CriTicaAL FLoopInG CONDITIONS 


Equation (4) can be expressed explicitly in G2 as 


(6) 


On differentiation this gives 


Hence, G? has a maximum or minimum value 
‘ 

when Bp? = 8, i.e. when ( °) = 2 which 

is equivalent to ee)" 16. 


(7) 


As can be seen 
4 
from Fig. 1 this value of - °) is well outside 
L 
the possible operating range of packed columns. 
Further differentiation of equation (7) leads to 
the following expression 


d B 


= — (9 — Bod) 


8 
64A (8) 
On substitution of Bp = 8 into the above 
equation it is found that 


(Gp) 
Since the value of the second differential 


coefficient is negative at (4) = 2, the 


function G2 = (pg) is then a maximum. This 
maximum is, however, outside the normal 
operating range so that, in practice, G? and thus 
G, always increases with increasing pg. Hence, in 
a packed column operated at reduced pressure, 
flooding would be expected to start at the top of 
the column where the pressure and therefore the 
density of the gas phase are lowest. The value 
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G3? = Beat, 
where B = 4 (4) 
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of the flooding rate to be used for design purposes of the vacuum pump and not at the pressure 
has, thus, to be calculated at the suction pressure existing in the still. 
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Letters to the Editors 


Discussion of the paper ‘‘ The partial coefficient of heat-transfer in 
a drying fluidized bed’’ by P. M. Heertjes and S. W. McKibbins* 


( Received 8 October 1956) 


Tue above paper concerns the measurement of heat-transfer 
coefficients in a fluidized bed of wet silica gel particles. 
Reference is made to similar work by KeErTTeNnRING, 
MANDERFIELD, and Smrru (Chem. Eng. Progress 1950 46 
139) and to heat transfer from air to coal without drying 
by Watton, Ovsen, and Levensprer (Ind. Eng. Chem. 
1952 44 1474). 

Apparently the authors were unaware of the work by 
the writer and colleague (WamsLEY and JoHanson, Chem. 
Eng. Progress 1954 50 347) in which heat-transfer coeffi- 
cients from gas to fluidized solid were obtained using an 
unsteady-state technique, and in which the results of 
KETTENRING, MANDERFIELD, and Surrn were extensively 
discussed. 

The heat-transfer coefficients obtained by the unsteady- 
state technique were only about one-eighth as large as 
those of Kerrenrinc, MANDERFIELD and Sairu, or WaAL- 
Ton, O_sen, and Levenspre.. The data of Ketrenrine, 
MANDERFIELD, and Smirn were reconciled with the 
unsteady-state data by assuming the solid temperature in 
the drying experiments was equal to the adiabatic satura- 
tion temperature: of the incoming air, rather than the 
temperature used by the authors. 

Heertses and McKissins, using an experimental 
technique developed earlier by Heertses, Boer, and 
DE Haas VAN Dorsser (Chem. Eng. Sci. 1953 2 97) were 
able to obtain results seemingly more reliable than those 
of Kerrenrinc, MANDERFIELD, and Smrru. The transfer 
coefficients resulting were stated by the authors (HEERTJES 
and McK ussrns) to differ from those of Kerrenrino, et al., 
and to agree more nearly with those of WaLTon ef al. All 
three sets of transfer coefficients are of the same magnitude, 
however, and therefore are about eight times as large as 
the values found by the unsteady-state technique of 
Wams .ey. In view of the increasing industrial applications 
of the fluidized-bed technique, such a large discrepancy in 
transfer coefficients must be reconciled, either by further 
work, or by showing the interpretation of data by one or 
the other technique to be faulty. 

In both researches involving wetted particles, the solid 
temperature during fluidization was not known. It was 
assumed equal to the exit gas temperature in each work. 
The following questions come to mind regarding this 
assumption : 

(a) Why would not the solid temperature approach 

adiabatic saturation (or wet bulb) temperature in 
such an adiabatic drying bed ? In Fig. 2 of HEERTJEs 


* Chem. Eng. Sci, 1956 5 161 


and McKissins, the temperature of the gel and of 
exit air is labelled as 15-4°C. The saturation 
temperature for the inlet air used is about 13°F, 
or — 10°C. 

If solid and gas reach equilibrium temperature 
within a few millimetres of the bed entrance, why 
would not the exit air be saturated with moisture 
at this temperature? In Fig. 11 of HEeerrtses, 
pe Boer, and pe Haas van Dorsser, the exit gas 
is shown to have a dewpoint of 29-8°F. If this air 
were in equilibrium with wet solid at a temperature 
of 12-1°C, as indicated by the authors, its saturation 
temperature and dewpoint should also be 12-1°C 
or 53-8°F. 

These considerations would seem to indicate that the 
leaving air is not in equilibrium with the wet solid particles. 
Referring to Fig. 10 of the 1954 paper of Wamsiey and 
JOHANSON, it will be noted that Wams_ey was able to 
find a finite difference between solid temperature and exit 
gas temperature, using the same technique for measuring 
solid temperature employed by Heertses and McK 
At the conditions of the latter work, this difference would 
be smaller and more difficult to detect. 

It is of importance to note that the major question is 
not the magnitude of this difference, but rather whether 
such a difference in solid and gas temperature exists 
throughout the bed. The existence of even a small driving 
force throughout the bed, rather than through only the 
first six millimetres as in Fig. 2 of Heertses and 
McK ssrns paper, has a profound effect on the magnitude 
of the transfer coefficient. Assuming a total bed depth of 
60 mm, consistent with Herertses’s earlier work, the 
transfer coefficient resulting from the data of Fig. 2 would 
be decreased about threefold if the solid temperature were 
15°C, rather than 15-4°C as assumed. 

The above Fig. 2 indicates a large gradient of gas 
temperature in the bottom of the bed, approaching exit 
conditions only 6 mm above the bed support. The writer 
interprets this gradient not as an indication of the rate 
of heat and mass transfer, but simply as a manifestation 
of the mixing of inlet gas with the gas already in the 
bed. The lack of approach to equilibrium of the gas and 
solid can be attributed to the bypassing of large bubbles 
or pockets of gas, which do not contact the solid effectively. 


University of Washington 
Seattle, Washington, U.S.A. 


L. N. JoHaNsoNn 
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Answer by P. M. Heertjes 


A POINT-WISE answer to the discussion by L. N. Jonanson 
of the paper: The Partial Coefficient of Heat Transfer in 
a Drying Fluidized Bed (Chem. Eng. Science 1956 5 161) 
seems indicated. 

(1) The authors were indeed unaware of the work of 
Wamsiey and Jonanson (Chem. Eng. Progress 1954 50 
847), which they regret. 

(2) Whether air leaving a wet fluidized mass can reach 
the adiabatic saturation temperature of the incoming air 
will depend largely on the water content of the gel. If 
the water content is such that the water at the surface 
of the gel behaves as free water and if the apparatus is 
really adiabatic, the adiabatic saturation temperature of 
the incoming air will be reached if the bed is high enough, 
viz. in the experiments presented higher than about 1 em. 

In the experiments under discussion however in the 
“ constant rate " period the water content of the gel was 
such that the sorbed water does not behave as free water 
and adsorption phenomena and energies have a strong 
bearing on the temperature of the leaving air. 

It can be found from the data presented that the 
humidity of the leaving air does not differ more than the 
experimental error, from the equilibrium humidity of the 
gel at the temperature of the leaving air (which in our 
opinion under the conditions given is equal to the gel- 
temperature). In general it is not correct to use calculated 
values based on the behaviour of free water as the tem- 
perature of the gel and it is incorrect in the cases 
evaluated. 

(3) The equations as used by Wamsiey and Jonanson 
will apply for any zone of heat interchange if for A will 
be taken the real surface of that zone. Wamsiey and 
JOHANSON take for A the surface of the whole bed, without 
giving real proof that this is allowed, e.g., by measuring 
the air temperatures over the bed. Therefore, the measure- 
ments as carried out only give certainty as to the value 
of the product h,A. If the actual exchanging surface is 
taken eight times too large, h, will be found eight times 
too small. In view of the amount of solid taken, this is 
quite possible. 


(4) As the only proof that the whole surface of the 
fluidized mass is involved in the heat exchange, WamsLEY 
and JoHANSON have measured the temperature of the 
solid mass, by using the same technique, viz., shutting 
off the air temporarily, as used by the present writer 
et al. As has been mentioned by the present writer these 
measurements indicate an equality of gel and leaving air 
temperature. However it has also been said that the 
place of such a measurement and moreover the evaluation 
of the time-temperature curve is very important to decide 
whether any point of this curve will be near the tempera- 
ture of the solid mass. In the measurements of WamMsLEY 
and JOHANSON with non-steady state conditions there can 
also be a time lag. Moreover, in the cases presented by 
the present writer ef al., other indications have been 
obtained which show that heat exchange and also mass 
exchange (not yet published) only occurs over a relatively 
small zone. These are : equal temperature of the air over 
the greatest part of the bed, the constant humidity over 
the same zone, the equality of a mean temperature 
measured in the bed with comparison of the temperature 
of the air in this zone. 

Some backmixing may occur but the experiments with 
different amounts of gel, indicate (see Chem. Eng. Science 
1953 2 97) that this influence cannot be large. 

(5) In the work of Jonanson et al., many of the 
magnitudes used have been calculated based on certain 
assumptions. In the Delft work material and heai balances 
could be drawn and checked with measured magnitudes. 
Including the heat of adsorption (to be found in the 
literature) besides the heat of evaporation the balances 
tallied within 5 to 10%. 

(6) It may finally be remarked that the coefficient as 
calculated by the present writer certainly does not repre- 
sent a correct value, because it has been calculated on a 
mean porosity over the whole bed. Work is carried out 
om now. P. M. Heerrses 
Laboratory for Chemical Engineering of the Delft University 
Delft 
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Book reviews 


P. ©. Carman: Flow of Gases through Porous 
Media. Butterworths, 1956. 182 pp. 30s. 


Tue flow of gases through porous media is of importance 
in many operations of interest to the chemical engineer, 
and Dr. CarMan has written a useful monograph on this 
rapidly developing field. Viscous flow through uncon- 
solidated beds of “normal” or “abnormal” types is first 
considered in the light of the Kozeny equation, the latter 
type implying departure of flow behaviour from the 
equation owing to non-uniform pore texture and high 
tortuosity. The author then deals with departures from 
Poiseuille flow when the pore diameter is comparable with 
or less than the mean free path, and discusses slip and 
Knudsen flow : with sufficiently fine pores such types of 
flow are possible at atmospheric pressures. This theme is 
later developed to include the effects of surface diffusion 
along the capillary walls ; this promises to be a method 
for isotopic separation. A chapter is devoted to the 
measurement of surface area by permeability studies, 
with a very fair discussion of the rival merits of nitrogen 
adsorption and permeability for this purpose. It should 
be remembered that for many materials, such as cement 
or textile fibres, the nitrogen method is less valuable 
since it gives the total surface, including the internal. 


There is a short section on gas separation, including 
chromatography (which most readers would prefer to study 
in separate monographs), turbulent flow and fluidisation 
of beds. 


The treatment is rather condensed and perhaps could 
be improved by a more detailed reference to experimental 
techniques. It is unfortunate that much of the theory is 
inevitably semi-empirical, which at times gives the reader 
a feeling of intellectual discomfort. This might have been 
relieved by a more fundamental discussion of those parts 
for which a rigorous theory exists, e.g. Knudsen flow. 
There are some omissions, e.g. Becker's striking work 
on surface diffusion studies by the use of the field emission 
microscope. 


These, however, are minor criticisms, which probably 
stem from the author's desire to keep down the size of the 
monograph. Dr. CARMAN’s writes with authority in a field 
to which he has contributed many papers, and is to be 
congratulated on bringing together much recent work on 
such diverse fields as soil mechanics and molecular sieves. 
One pleasing feature of a technical book is the use through- 
out of ¢.g.s units. The production is good and the price 


is reasonable. R. S. Brap.ey 


R. S. Antes and R. D. Newron : Chemical Engineer- 
ing Cost Estimation. McGraw-Hill Rook Company, 
1955, pp. 263, $6.00 (45s.). 


To present a simple system for arriving at a cost estimate 
for chemical plant, of sufficient accuracy to enable mana- 
gerial decisions to be made, seems an impossible aim. This 
book attempts just that. It is a revised version of the 
book of the same name in the McGraw-Hill Chemical 
Engineering Series. No doubt a skilled cost accountant 
might point out flaws in this book and no doubt the many 
plots of cost versus capacity for a wide range of unit 
operations equipment, piping, insulation and so on, 
represent a gross over simplication. Nevertheless, to the 
chemical engineer, particularly in the smaller firm, required 


to produce quickly an estimate of the cost of a projected 
operation, this book can be recommended as a mine of 
useful information, clearly presented. 

The book is entirely centred on the U.S. market. This 
means that if a direct conversion of U.S. dollars to local 
currency is made at a convenient rate of exchange, many 
of the costs so found will be useless outside of the United 
States. It is a serious criticism of the book, from the 
point of view of the European reader, that this problem 
is completely ignored, and no solution is hinted at. An 
adidtional chapter on means suggested for convering costs 
to other currencies should be seriously considered if the 
book is to be of real use outside of the U.S.A. 


D. M. Wits0on 
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